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I.  INTRODUCTION 

A  growing  number  of  defense  applications  of  microwave  and  millimeter 
wave  radiation  present  requirements  that  cannot  be  satisfactorily  met  with 
existing  fundamental  mode  sources.  These  requirements  include  high  power 
at  high  frequency  with  high  efficiency,  longer  pulse  lengths,  tunability,  and 
amplifier  operation  together  with  compact  size,  ruggedness  and  good 
reliability.  An  illustrative  example  in  electronic  warfare  is  the  destruction 
of  electronic  controls,  sensors,  and  communication  systems.  There  the  most 
effective  approach  to  coupling  microwave  power  into  the  target  utilizes 
super-power  pulses  with  rapid  rise  times  at  the  characteristic  operating 
frequency  of  the  target  system.  These  frequencies  are  already  high  and 
expected  to  be  higher  in  next-generation  systems.  Even  now,  high  power  is 
available  only  at  frequencies  below  the  characteristic  frequencies  of 
offensive  weapons  systems;  and  coupling  must  therefore  be  primarily 
through  apertures  not  intended  for  microwave  propagation  with  far  lower 
efficiency  than  would  be  possible  for  direct  coupling. 

The  capabilities  of  existing  sources  for  meeting  these  requirements  are 
limited  largely  because  the  generation  of  power  is  based  on  the  passage  of 
high-energy  electron  beams  through  resonant  structures.  The  maximum 
power  emitted  by  such  devices  decreases  rapidly  with  increasing  frequency, 
Pf2  ~  constant,  where  P  is  the  power  emitted  at  frequency  f.  This  limitation 
arises  principally  because  the  size  or  the  resonant  structure  Is  proportional 
to  the  wavelength,  A  =  c/f;  and  increasing  the  frequency  (decreasing  the 
wavelength)  increases  the  power  density  that  must  be  handled  by  the 
resonant  structure.  The  resulting  power  density  limitation  in  conventional 
microwave  sources  is  further  exacerbated  at  high  power  by  the  need  to 
increase  the  electron  beam  voltage  to  reduce  space-charge  effects  which 
would  otherwise  degrade  the  source  performance. 
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The  microwave  source  concept  proposed  here  circumvents  the 
conventional  beam  transit  time,  resonant  structure,  and  space  charge 
limitations  by  storing  energy  in  a  magnetic-mirror  confined,  electrically 
neutral,  hot-electron  plasma  and  periodically  transforming  a  fraction  of  this 
stored  energy  into  pulses  of  very  high  peak  power.  The  transformation  is 
mediated  by  a  whistler  wave  launched  at  one  end  of  the  mirror-confined 
plasma,  propagating  along  the  magnetic  lines  of  force,  amplified  by  the 
anisotropic  hot-electron  plasma,  and  subsequently  collected  and  focused  by 
quasi-optical  elements  at  the  opposite  end  of  the  plasma.  The  operating 
frequency  of  this  amplifier  is  determined  by  the  strength  of  the  static 
magnetic  field  and  can  therefore  be  scaled  to  very  high  values.  The  gain  is 
determined  by  the  physical  dimensions  and  the  energy  density  of  the 
mirror-confined  plasma  so  that  the  output  power  can  be  increased  by  scaling 
in  size.  Because  both  the  operating  frequency  and  the  plasma  energy  density 
increase  with  increasing  magnetic  field  strength,  it  appears  that  the  peak 
power  achieved  in  the  proposed  concept  will  increase  with  frequency  (as 
well  as  the  volume  of  the  magnetic  mirror)  to  give  highly  favorable  scaling 
properties  for  defense  applications.  The  resulting  source  can  operate  in 
amplifier  or  oscillator  modes,  can  have  high  efficiency  and  be  implemented 
in  a  rugged,  compact  embodiment. 

The  mirror-confined,  hot-electron  plasma  that  comprises  the  active 
medium  responsible  for  amplification  of  the  incident  whistler  results  from 
the  collective  interaction  of  two  or  more  oppositely  directed  electron 
streams  with  a  mirror-confined  background  plasma.  Experimental  and 
theoretical  studies  of  hot-electron  plasmas  created  by  beam-plasma 
interactions  confirm  the  presence  in  the  plasma  of  many  separate  regions  of 
high  field  strength  which  are  individually  coherent  but  mutually  incoherent. 
Such  a  volume  distribution  of  independent  regions  is  not  suitable  for  the 
direct  generation  of  useful  microwave  power;  but  it  is  extremely  efficient 
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for  stochastic  heating  of  a  minority  of  the  plasma  electrons  to  create  a 
mirror-confined  hot-electron  plasma.  The  energy  stored  in  this  plasma  can, 
in  turn,  be  used  either  to  amplify  incident  whistler  waves  propagating  along 
the  lines  of  force  or  as  the  source  of  spontaneous  oscillations  of  whistler 
wave  power.  In  effect,  energy  is  stored  throughout  the  plasma  volume  by  the 
beam-plasma  interaction  but  transformed  into  a  coherent,  two-dimensional 
form  by  the  unstable  whistler  wave. 

To  recapitulate,  the  high-power  microwave  source  concept  to  be 
developed  in  this  project  circumvents  many  of  the  limitations  inherent  in 
conventional  microwave  sources: 

1)  Limitations  arising  from  the  finite  transit  time  of  beam  electrons 
through  the  resonant  structure  are  circumvented  by  using  the 
energy  stored  in  a  relativistic-electron  plasma  that  is  confined  in 
a  magnetic-mirror  trap. 

2)  Space-charge  limits  on  the  beam  are  removed,  since  the 
relativistic  electron  plasma  is  electrically  neutral. 

3)  The  need  for  a  metallic  slow  wave  structure  is  obviated  by  having 
the  stored  energy  transformed  into  microwave  fields  through  the 
amplification  of  whistler  waves  that  are  guided  by  the  magnetic 
field  and  the  plasma. 

4)  The  technology  for  storing  energy  in  the  hot-electron  plasma 
uses  simple  steady-state  streams  of  electrons  flowing  along  the 
magnetic  lines  of  force  through  a  background  plasma,  and  can  be 
implemented  in  a  rugged,  compact  embodiment. 

5)  The  operating  frequency  of  the  source  is  determined  by  the 
strength  of  the  static  magnetic  field,  and  can  therefore  readily 
be  tuned  as  well  as  scaled  to  high  values. 

6)  The  gain  and  maximum  output  power  are  determined  by  the 
physical  dimensions  and  the  energy  density  of  the  mirror-confined 
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hot-electron  plasma,  so  that  output  power  can  be  increased  by 
scaling  in  size,  and  external  magnetic  field  strength. 

Results  achieved  in  Phase  I  of  this  project  and  summarized  in  Section 
3.1  confirm  the  feasibility  of  accumulating  significant  quantities  of  kinetic 
energy  in  a  mirror-confined  hot-electron  plasma  using  simple 
interpenetrating  streams  of  electrons  in  a  background  plasma,  as  indicated 
schematically  in  Figure  1.  Under  suitable  operating  conditions, 
electrostatic  waves  propagating  nearly  parallel  to  the  magnetic  lines  of 
force  are  driven  to  large  amplitude  by  the  collective  interaction  of  the 
electron  streams  with  the  cold,  background  plasma.  (1)  Nonlinear  processes 
are  known  to  couple  the  large  amplitude  electrostatic  waves  to 
electromagnetic  waves  at  twice  the  local  upper  hybrid  frequency.  (2)  These 
waves  propagate  outward  into  regions  where  they  can  be  preferentially 
absorbed  by  the  energetic,  mirror-confined  electrons  through  cyclotron 
damping  at  overtones  of  the  relativistic  electron  gyrofrequency.  (3)  If  the 
plasma  is  enclosed  in  a  suitable  conducting  shell,  the  electromagnetic  waves 
can  be  confined  in  the  plasma  region  and  absorbed  by  the  hot-electron  plasma 
in  ways  that  are  fundamentally  identical  to  the  processes  responsible  for  hot 
electron  generation  in  microwave  heated  plasmas.  Other  electrostatic  waves 
propagating  at  oblique  angles  to  the  magnetic  field  can  also  be  driven  to 
large  amplitude  by  the  beam-plasma  interaction.  These  waves  extract 
parallel  kinetic  energy  from  the  electron  streams  and  increase  their 
perpendicular  kinetic  energies.  As  a  result,  some  of  the  energetic-beam 
electrons  are  trapped  in  the  magnetic  mirror. 

The  mirror-confined  hot-electron  plasma  resulting  from  the 
beam-plasma  interaction  forms  an  active  medium  that  can  be  used  to  amplify 
incident  whistler  waves  and  thereby  transform  the  kinetic  energy  of  the 
hot-electron  plasma  into  microwave  field  energy.  An  indication  of  the 
plasma  conditions  under  which  high-gain  amplification  of  whistlers  can 
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occur  is  shown  in  Figure  2,  where  the  threshold  for  growth  is  plotted  as  a 
function  of  the  hot-electron  temperature  anisotropy  and  relative  energy 
density  (4),  as  discussed  at  length  in  Section  3.3.  For  efficient  operation  of 
this  concept  it  is  necessary  to  accumulate  the  hot-electron  plasma  under 
stable  conditions,  and  then  rapidly  alter  the  plasma  state  so  that  whistler 
amplification  is  possible.  This  rapid  alteration  in  the  plasma  conditions  can 
be  accomplished  in  principle  through  adiabatic  magnetic  compression, 
reduction  of  the  cold-plasma  density,  or  a  combination  of  the  two.  Adiabatic 
magnetic  compression  appears  to  be  particularly  well  suited  to  amplifier 
operation,  while  supplemental  electron  cyclotron  heating  for  control  of  the 
hot-electron  temperature  anisotropy  and  the  cold-plasma  density,  together 
with  the  use  of  quasi-optical  resonator  techniques,  is  especially  attractive 
for  operation  in  the  pulsed  oscillator  mode.  Both  the  oscillator  and  amplifier 
modes  are  expected  to  yield  high-power  pulses  of  coherent  radiation.  (5) 

1.1  Identification  of  the  Key  Issues 

Earlier  work  (1)  on  the  formation  of  hot-electron  plasmas  by 
beam-plasma  interactions  has  not  included  an  adequate  characterization  of 
the  equilibrium  properties  of  the  hot-electron  component.  In  order  to  utilize 
these  hot-electron  plasmas  to  amplify  whistler  waves  it  is  necessary  to 
characterize  and  control  the  spatial  distribution  of  hot-electrons  as  well  as 
their  temperature  anisotropy  and  relative  energy  density.  Techniques  for 
determining  the  most  important  hot-electron  plasma  parameters  have  been 
developed  in  connection  with  electron  cyclotron  heating  research,  (6)  and 
these  techniques  can  be  applied  to  the  hot-electron  plasma  generated  by  the 
beam-plasma  interactions. 

The  optimization  of  beam-plasma  generated  hot-electron  plasmas  for 
the  PEMS  application  can  best  be  determined  through  an  experimental 
demonstration  of  whistler  wave  gain,  for  operation  in  the  amplifier  mode,  or 
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triggered  high-power  microwave  pulses,  for  operation  in  the  oscillator  mode. 

With  regard  to  hot-electron  plasma  formation  and  whistler 
amplification,  a  valid  theoretical  model  is  needed  to  facilitate 
interpretation  of  experiments  and  optimization  of  the  concept  as  well  as 
scaling  the  results  of  Phase  II  experiments  to  Phase  III  conceptual  devices. 
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(1]  DC  Power  input  via  steady-state  electron  beams 

(2]  Electromagnetic  power  at  2  go,*  (-20^)  generated  by  beam-plasma 
interactions 

(3)  Power  generated  by  beam-plasma  interactions  stored  in  anisotropic, 
mirror-confined,  hot-electron  plasma 

(4)  Coherent,  pulsed,  high-power  microwave  output  from  unstable 
whistler  waves 


FIGURE  la:  Schematic  representation  of  the  Two-Stream  PEM5  Concept 
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growth  of  whistler  waves  in  spatially  inhomogeneous, 
anisotropic,  finite-beta  hot-electron  plasmas 
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2.  TECHNICAL  APPROACH 


2.1  Technical  Objectives 

In  order  to  demonstrate  experimentally  the  potential  advantages  of  the 
Two-Stream  PEMS  concept  enumerated  in  the  preceeding  section,  we  have 
developed  an  experimental  approach  that  can  be  implemented  in  a  timely  and 
cost-effective  way  using  the  existing  AMPHED  facility  at  AMPC,  described  in 
Section  3.2.  This  approach  seeks  to  achieve  three  major  objectives: 

1 )  an  experimental  demonstration  and  characterization  of  efficient 
energy  accumulation  using  the  two-stream  beam-plasma 
interaction  in  suitable  conducting  cavities; 

2)  the  experimental  demonstration  of  whistler  amplification  by 
the  mirror-confined  hot-electron  plasma  and/or  triggered, 
high-power  pulses  in  the  oscillator  mode  of  operation;  and 

3)  the  development  of  theoretical  models  to  interpret  the 
experimental  results  and  evaluate  the  performance  of 
conceptual  systems  using  the  Two-Stream  PEMS  approach. 

These  objectives  can  be  accomplished  through  an  experimental  program 
that  utilizes  the  existing  AMPHED  experimental  facility  together  with  a 
closely  coupled  theoretical  effort  aimed  at  developing  a  valid  theoretical 
model  of  the  experimental  results.  In  this  section  we  describe  in  some 
detail  the  technical  approach  to  be  used  in  this  project. 
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2.2  Modification  of  the  AMPHED  Facility 

The  existing  AMPHED  experimental  facility  requires  only  minor 
modifications  to  permit  the  experimental  production  and  characterization  of 
hot-electron  plasma  generated  by  beam-plasma  interactions.  A  symmetrical 
structure  to  generate  two  oppositely  directed  electron  streams  will  be 
located  in  AMPHED  as  shown  in  Figure  lb.  The  streams  will  initially  be 
powered  by  an  existing  10  kV,  500mA  power  supply.  The  flat-field  magnetic 
configuration  in  AMPHED  provides  a  uniform  interaction  region  up  to  80  cm 
long.  This  large  region  of  uniform  magnetic  intensity  can  be  used  to  ensure 
strong  growth  of  unstable  electrostatic  waves.  The  background  plasma  can 
be  generated  spontaneously  by  beam  ionization  of  the  ambient  gas,  in  which 
case  the  major  control  of  background  plasma  density  is  through  control  of 
the  ambient  gas  pressure.  In  addition,  since  AMPHED  is  supplied  with  several 
sources  of  microwave  power  suitable  for  creation  and  control  of  the 
background  plasma  density  and  temperature,  these  can  be  used  to  achieve  an 
unusual  degree  of  control  over  these  parameters  for  studies  aimed  at 
optimum  generation  of  PEMS-relevant  hot-electron  plasmas. 

Control  of  the  background  plasma  density  is  expected  to  be  particularly 
important  in  optimizing  hot-electron  production,  since  the  electromagnetic 
power  generated  by  the  non-linear  electrostatic  waves  is  mainly  at 
frequencies  near  twice  the  local  upper  hybrid  frequency.  We  anticipate  the 
most  efficient  generation  of  hot  electrons  if  this  power  is  near  the  second 
harmonic  of  the  relativistic  electron  gyrofrequency  in  the  region  where  hot 
electrons  accumulate.  Control  of  the  background  plasma  density  will  allow 
us  to  vary  the  upper  hybrid  frequency  to  test  this  hypothesis. 

The  DC  electrical  supply  for  the  AMPHED  magnetic  coils  has  recently 
been  upgraded  to  provide  the  level  of  control  required  for  precise 
determination  of  the  diamagnetic  fields  created  by  hot-electron  plasmas 
using  arrays  of  Hall  probes.  (6)  The  spatial  dependence  of  the  hot-electron 
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energy  density  will  be  measured  experimentally  with  these  Hall  probe  arrays 
as  well  as  by  the  traditional  skimmer  probes  used  successfully  to 
characterize  high  energy  density  hot-electron  plasmas  in  ELMO.  (3)  The  local 
distribution  of  hot  electrons  in  energy  will  be  measured  using  x-ray 
techniques  as  in  numerous  earlier  ECH  hot-electron  experiments.  Background 
plasma  parameters  will  be  determined  with  conventional  Langmuir  probes 
presently  in  use  on  AMPHED. 

2.3  Experimental  Characterization  of  Hot-Electron  Plasma 

Generation  Using  Beam-Plasma  Processes 

Following  the  installation  of  the  beam-forming  structure  we  anticipate 
an  initial  series  of  experiments  aimed  at  optimizing  the  beam-plasma 
generation  of  hot-electron  plasmas  and  establishing  the  operating  conditions 
for  which  satisfactory  hot-electron  production  is  achieved.  The  system 
parameters  expected  to  be  particularly  important  to  hot-electron  production 
are  the  electron  beam  energy  and  density,  the  background  plasma  density,  and 
the  spatial  extent  and  uniformity  of  the  two-stream  interaction  region.  We 
will  optimize  the  hot-electron  diamagnetic  energy  density  with  respect  to 
these  control  parameters  and  seek  to  determine  preliminary  scaling 
relations  that  can  be  used  to  extrapolate  the  Phase  II  results  for 
conceptualization  of  Phase  III  applications. 

In  order  to  facilitate  the  development  of  valid  theoretical  models  of  the 
hot-electron  plasma  formation,  we  will  monitor  the  power  spectrum  of  the 
electromagnetic  radiation  generated  by  the  large  amplitude  electrostatic 
waves  driven  by  the  beam-plasma  interaction.  In  this  way  we  will  seek  to 
estimate  the  overall  efficiency  of  the  process  as  well  as  determine  ways  of 
enhancing  that  efficiency. 
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Having  delineated  a  satisfactory  range  of  operating  conditions  for 
generating  hot-electron  plasmas,  we  will  undertake  a  more  detailed 
characterization  of  these  hot-electron  plasmas.  This  characterization  will 
include  the  spatial  properties  of  the  hot-electron  component  as  well  as  the 
local  distribution  of  hot  electrons  in  energy.  Conversion  of  the  hot-electron 
energy  into  microwave  power  depends  on  the  hot-electron  temperature 
anisotropy  and  relative  energy  density.  We  will  seek  to  identify  the 
dependence  of  these  properties  on  the  externally  controlled  operating 
parameters  such  as  ambient  gas  pressure  and  electron  stream  power,  as  well 
as  supplemental  ECH  and/or  UORH  power. 

2.4  Measurements  of  Whistler  Amplification 

At  the  present  time  AMPC  is  developing  experimental  techniques  to 
measure  the  propagation  and  growth  of  whistler  waves  launched  along 
magnetic  lines  of  force  and  passing  through  magnetic-mirror  confined 
plasmas.  Whistlers  are  launched  from  a  dielectric-loaded  wave  guide 
structure  whose  radiation  pattern  is  modified  by  a  shaped  dielectric  insert, 
as  described  in  Section  3.2.  Wavelengths  of  the  propagating  whistler  are 
measured  by  microwave  interferometric  techniques,  and  local  intensities  are 
determined  by  RF  probes  as  described  in  Section  3.2.  In  the  second  phase  of 
the  experimental  program  proposed  here,  we  will  utilize  similar  techniques 
to  help  characterize  the  amplification  of  whistlers  by  the  hot-electron 
plasmas  generated  through  the  beam-plasma  interactions.  The  amplified 
whistler  will  be  detected  and  its  power  level  determined  using  an  array  of 
horns  mounted  on  the  end  opposite  the  launch  structure. 
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2.5  Development  of  Valid  Theoretical  Models 

In  order  to  Interpret  the  experimental  results  to  be  obtained  in  Phase  II 
as  well  as  to  evaluate  conceptual  systems  to  be  considered  in  Phase  III, 
AMPC  will  develop  and  validate  a  theoretical  model  that  encompasses  the 
two  major  aspects  of  the  Two-Stream  PEMS  concept;  namely, 

1)  Efficient  Formation  of  Mirror-Confined  Hot-Electron  Plasmas  by 
Beam-Plasma  Processes;  and 

2)  Transformation  of  Plasma  Energy  to  Coherent  Microwave  Power 
Through  Whistler  Amplification. 

Within  the  first  category  we  will  seek  to  model  the  dominant  electrostatic 
beam-plasma  interactions,  the  processes  that  saturate  the  growth  of  the 
electrostatic  waves,  the  three-wave  processes  governing  the  generation  of 
electromagnetic  waves,  and  the  preferential  absorption  of  these  waves  by 
mirror-confined  hot  electrons.  The  basic  linear  dispersion  relation  code 
describing  growth  of  electrostatic  beam-plasma  modes  has  been  used 
extensively  in  Phase  I,  as  described  in  Section  3.1.  The  theory  of  the 
three-wave  coupling  has  been  developed  by  several  authors  (2)  and  can  be 
validated  using  Phase  II  experiments.  The  preferential  absorption  of 
high-frequency  electromagnetic  waves  by  mirror-confined  hot  electrons  has 
been  studied  experimentally  and  theoretically  by  AMPC  scientists  in  the  past 
few  years.  (3)  These  three  separate  components  will  be  linked  together  to 
form  the  principal  framework  for  the  first  major  aspect  of  the  Phase  II 
modelling  program. 

The  second  major  aspect  will  build  on  the  earlier  work  done  at  AMPC  to 
evaluate  realistic  threshold  conditions  and  maximum  amplification  rates  for 
whistler  waves  propagating  through  anisotropic,  mirror-confined, 
hot-electron  plasmas.  (4)  New  theoretical  studies  will  be  undertaken  to 
clarify  the  roles  of  relevant  non-linear  mechanisms  controlling  the  onset  of 
whistler  growth,  such  as  the  relative  cold-plasma  density  and  temperature. 
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The  mechanisms  responsible  for  saturating  the  growth  of  whistlers  will  also 
be  investigated  in  order  to  estimate  the  maximum  output  power  in  the 
amplifier  mode. 

2.6  Summary  of  the  Phase  I  Technical  Objectives 

The  main  objective  of  the  Phase  I  program  was  a  preliminary  design  of 
an  experimental  device  capable  of  demonstrating  the  validity  of  the  basic 
elements  of  the  concept: 

1.  Effective  accumulation  of  significant  levels  of  energy  density  in  a 
mirror-confined,  hot-electron  plasma  using  interpenetrating 
electron  streams  in  a  background  plasma;  and 

2.  Amplification  of  externally  launched  whistler  waves  propagating 
along  magnetic  lines  of  force  through  the  mirror-confined 
hot-electron  plasma. 

In  order  to  accomplish  this  design  objective,  it  was  necessary  to 
address  a  number  of  fundamental  issues  whose  clarification  provided  a  solid 
theoretical  framework  for  the  proposed  concept.  Among  the  most  immediate 
issues  were  the  following: 

1 .  Optimal  Conditions  for  Creating  the  Hot-Electron  Plasma 

A  large  empirical  base  of  information,  together  with  a  body  of 
somewhat  idealized  theory  presently  exists  regarding  the  use  of 
beam-plasma  interactions  to  create  mirror-confined  hot-electron  plasmas. 
This  material  was  updated  using  relevant  results  from  experiments  and 
Fokker-Planck  models  of  the  electron  cyclotron  heating  (ECH)  process  which 
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clearly  establish  the  role  of  O-mode  fields  at  overtones  of  the  hot-electron 
gyrofrequency  in  forming  high-beta,  hot-electron  plasmas.  In  this  way,  we 
sought  to  determine  the  available  range  of  operating  parameters  and  the 
sensitivity  of  the  electron  heating  process  to  the  major  device  parameters 
entering  the  design. 

2.  Alternative  Ways  of  Controlling  the  Release  of  Stored  Energy  for 
High-Gain  Amplification. 

Two  general  approaches  to  bringing  the  plasma  rapidly  to  the  threshold 
for  whistler  amplification  are  presently  identified;  namely,  control  of  the 
hot-electron  temperature  anisotropy  and  the  magnetic  field  uniformity,  and 
control  of  the  background  cold-plasma  density  and  temperature.  A  more 
detailed  analysis  of  each  approach  is  needed  to  implement  one  or  both  of 
these  approaches  in  the  device  design. 

3.  Quasi-Optical  Output  Coupling  Structures 

An  analysis  of  the  propagation,  growth  and/or  damping  of  whistler 
waves  in  candidate  device  concepts  is  needed  to  support  the  design  of 
effective  coupling  structures  to  guide  the  amplified  whistler  wave  into 
suitable  microwave  transmission  systems. 


3.  MAJOR  RESULTS  ACCOMPLISHED 

3.1  Creation  of  Hot-Electron  Plasmas  Using  Beam-Plasma 
Interactions 

In  Phase  I  of  this  project  we  have  investigated  the  fundamental 
processes  expected  to  play  dominant  roles  in  determining  the  efficiency  with 
which  beam-plasma  interactions  will  generate  a  substantial  density  of 
relativistic,  mirror-confined  electrons.  These  include  three  major 
phenomena; 

1)  the  growth  of  unstable  electrostatic  waves  driven  by  the  two 
streams  of  electrons  interacting  with  the  background  plasma; 

2)  the  non-linear  coupling  of  oppositely-directed  pairs  of 
large-amplitude  electrostatic  waves  to  excite  an 
electromagnetic  wave  propagating  radially  outward;  and 

3)  the  preferential  absorption  of  this  electromagnetic  wave  by 
energetic  electrons  confined  in  a  magnetic  mirror. 

Each  of  these  three  separate  phenomena  has  a  substantial  empirical  and 
theoretical  base  which  we  will  briefly  review  here.  Their  integration  to 
form  an  efficient  technique  for  accumulating  energy  in  a  mirror-confined 
hot-electron  plasma  is  a  major  objective  of  Phase  II  of  the  present  project. 
It  will  be  clear  from  the  review  of  each  separate  phenomenon  that  the 
overall  feasibility  of  the  concept  is  strongly  supported  and  the  prognosis  for 
efficient  energy  storage  by  the  proposed  technique  is  excellent. 
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3.1.1  Unstable  Electrostatic  Waves  in  the  Proposed 
Configuration 

The  proposed  proof-of-principle  plasma  configuration  consists  of  two 
co-linear  interpenetrating,  oppositely-directed  electron  streams  flowing 
along  the  axis  of  the  external  magnetic-coil  system  through  a 
low-temperature  background  plasma.  The  properties  of  unstable 
electrostatic  waves  in  this  configuration  have  been  studied  extensively  both 
experimentally  (I)  and  theoretically.  (7)  The  most  rudimentary  theoretical 
model  capable  of  describing  the  propagation  and  growth  of  electrostatic 
waves  in  this  configuration  consists  of  infinite,  homogeneous  plasma 
components  in  a  uniform  steady  magnetic  field.  The  dispersion  relation  for 
electrostatic  waves  in  such  a  plasma  was  derived  by  Harris  (8).  In  the  limit 
of  negligible  thermal  speeds,  this  dispersion  relation  was  shown  to 
correspond  well  with  laboratory  experiments  (I).  Here  we  summarize  the 
main  predictions  of  the  theoretical  model  bearing  on  the  excitiation  of  large 
amplitude  electrostatic  waves  that  can  interact  non-linearly  to  excite 
electromagnetic  waves  suitable  for  generating  a  mirror-confined 
hot-electron  plasma. 

The  dominant  mode  of  instability  arises  from  interactions  between  each 
separate  beam  and  the  background  plasma.  The  properties  of  the  dominant 
mode  are  thus  identical  to  those  of  the  single-beam  cold-plasma  system 
studied  by  O'Neil  et.  al.  (9)  The  frequency  of  the  fastest  growing  mode  is 
very  close  to  the  electron  plasma  frequency,  as  shown  In  Fig.  3.  Here  we 
display  the  wave  vector,  the  real  frequency,  and  the  growth  rate  of  the 
fastest  growing  mode  for  a  range  of  background  electron  plasma  frequencies 
with  the  beam  density  fixed.  The  wave  vector  has  been  chosen  to  maximize 
the  growth  rates.  Each  of  the  parameters  is  normalized  as  follows: 

1 )  the  wave  vector  is  expressed  through  the  dimensionless  quantity 
H  h  k,  ,vb/n,  where  k(|  is  the  component  of  the  wave  vector 
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parallel  to  the  static  magnetic  field,  vb  is  the  beam  speed,  and  fl 
is  the  electron  gyrofrequency; 

2)  the  wave  frequency  and  growth  rate  are  normalized  by  the 
electron  gyrofrequency,  UJ  =  a>/fi. 

3)  the  background  electron  plasma  frequency  is  normalized  by  the 
electron  gyrofrequency,  UJe  j  copg/n,  where  o^2  -  e?ne/m€0. 
The  electron  charge  and  mass  are  -  e  and  m,  the  background 
electron  density  is  and  €g  is  tho  permittivity  of  free  space, 

4)  the  electron  beam  density,  r>„  is  expressed  through  the 
equivalent  plasma  frequency,  IDt,  s  co^/fi,  with  co#2  =  e2r\,/m€o. 

As  is  clear  from  Figure  3,  the  fastest  growing  mode  satisfies  two 
approximate  criteria;  namely, 

Ui  We  or  co  cope 
and 

H  **  UJe  or  kvb  ^  cope- 

The  maximum  growth  rate  is  given  roughly  by 
UJj  ^  0.6tUJb2aje),/3, 

where  U>,  is  the  imaginary  part  of  the  frequency,  Y,  normalized  by  the 
electron  gyrofrequency.  Thus, 

Y  “  0.6(rvf\),/3cop,. 

Note  that  the  coefficient  0.6  is  close  to  the  corresponding  analytic  estimate 
(3,/2/Z4/3)  *  0.687  deduced  by  O’Neil  et.  al.  (9)  The  dependence  of  the  growth 
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Properties  of  the  fastest  growing  electrostatic  beam-plasma 
instability:  the  (parallel )  wave  vector  H=  kjj  v5/0;  the  real 
frequency,  UJr  =  co/Q;  and  the  growth  rate,  UJ,  =  y/O,  where  n 

is  the  electron  gyrofreuqency.  The  background  electron  plasma 
frequency  is 
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rate  on  the  electron-stream  density  is  shown  in  Fig.  4,  where  UJj  is 
displayed  as  a  function  of  UJb  for  UJe  =  1  and  H  =  1.04,  the  value  for  maximum 
growth  rate. 

In  obtaining  the  results  shown  in  Figs.  3  and  4  we  have  chosen  the  wave 
vector  to  be  parallel  to  the  magnetic  field,  i.e.,  kA  <*  0.  For  electron  streams 
of  finite  radial  extent  the  boundary  conditions  at  the  surface  of  the  stream 
may  require  a  finite  value  of  kA.  For  example,  a  suitable  boundary  condition 
for  many  cylindrical  systems  is  J0{kAa]  =  0,  where  J0  is  the  Bessel  funcion  of 
order  zero,  and  a  is  the  radius  of  the  column.  For  such  a  case, 
kA  *  Z4/a.  lo  illustrate  how  oblique  propagation  affects  the  growth  rate  of 
the  modes  shown  in  Fig.  3  and  4,  we  display  in  Fig.  5  a  plot  of  U1,[0],  where  0 
■  terrHkA/k||).  For  the  particular  case  chosen  here,  UJe=1,  11%=  0.1  and  H  = 
1.04,  0.8,  and  0.6.  The  growth  rate  is  substantially  reduced  for  0-0;  but  y 
vanishes  only  for  propagation  nearly  perpendicular  to  the  magnetic  field.  It 
is  clear  from  Fig.  5,  however,  that  the  obliquely  propagating  wave  is  entirely 
distinct  from  the  wave  that  propagates  parallel  to  the  magnetic  field. 
Indeed,  the  growth  rate  is  maximized  for  non-zero  angles  of  propagation  for 
H  =  0.8  and  0.6;  whereas  the  maximum  growth  rate  occurs  for  0  =  0  and  H  = 
1.04. 

The  sub-dominant  obliquely-propagating  mode  arises  from  an 
interaction  between  the  two  counter-streaming  electron  beams,  in  which  a 
longitudinal  space-charge  wave  on  one  stream  couples  to  a  cyclotron  wave  on 
the  opposite  stream.  A  particularly  noteworthy  feature  of  this  mode  is  that 
it  transfers  parallel  momentum  from  each  stream  into  perpendicular 
momentum  of  the  oppositely-directed  stream.  Excitation  of  this  mode  could 
significantly  enhance  the  magnetic  trapping  of  energetic-beam  electrons 
collectively  scattered  out  of  the  streams. 
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A.  Dependence  of  growth  rate  on  the  electron  beam  density 
UJb  =  t-Ope/ft  >s  the  electron  plasma  frequency  corresponding 
to  the  beam  density. 
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The  electrostatic  waves  are  predicted  to  grow  at  the  linear  growth  rate 
until  the  associated  electrostatic  potential  is  large  enough  to  trap  the 
be^m  electrons  (9)  The  growth  then  ceases  at  this  level;  namely, 


e^t  “  m(vb  -  o>/k)jF/2  =  fc^l  -  U^/H)2, 

where  »  mvb2/2  is  the  kinetic  energy  of  the  beam  electrons.  The 
parallel  phase  velocity  of  the  fastest  growing  wave  is  proportional  to  but 
slightly  less  than  the  beam  velocity.  In  Fig.  6  we  display  1-U!r/H  as  well  as 
other  wave  parameters  versus  (r\,/ne],/3)  from  which  we  obtain 

1  -  Uir/H  -  0.5(rv'ne),/3, 

in  agreement  with  O'Neil  et.  al.  (9)  When  the  linear  growth  saturates,  a 
fraction  of  the  kinetic  energy  of  the  beam  electrons  is  converted  to 
electrostatic  energy: 

t^/2  rx.Eblr^/rieF3/*. 

This  process  has  been  observed  in  numerous  simulations  of  the  beam-plasma 
interaction.  (10) 

O'Neil  et.  al,  (9)  show  how  the  evolution  of  the  fastest  growing  wave 
leads  to  a  very  narrow  spectral  distribution  of  the  wave  energy  in  wave 
number.  The  initial  width,  Akjj/k,,  ~  (r^/r^)1^3,  is  predicted  to  be  reduced 
to  a  much  narrower  range  given  by 

$k/k0  -  3{ln2/N)1/2  2_5/6(r\,/r^)1/3, 
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FIGURE  6;  Dependence  of  properties  of  the  fastest  growing  wave  on  the 
relative  density  of  the  beam  electrons  and  the  background 
plasma  electrons,  r>,/rv  Here  £  s  (r\,/i\),/3 
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where  N  is  the  number  of  e-foldings  in  wave  amplitude  before 
saturation  of  the  linear  growth,  N  =  OHO2).  The  spectral  density  of 
electrostatic  field  energy  is  predicted  to  increase  roughly  a  factor  of  7  -  10 
through  this  narrowing  of  the  spectral  distribution.  As  O'Neil  et.  al.  (9) 
observe,  the  saturated  wave  is  essentially  a  pure  sinusoidal  disturbance. 

3.1.2  The  Non-Linear  Three-Wave  interaction 

The  large  amplitude  electrostatic  waves  described  in  Sect.  3.1  are 
generated  by  each  of  the  co-1  inear,  oppositely-directed  electron  streams 
interacting  with  the  background  plasma.  When  pairs  of  these 
oppositely-directed  electrostatic  waves  interact  they  can  generate  an 
electromagnetic  wave  propagating  radially  outward  from  the  electron  beam 
region.  This  fusion  of  two  electrostatic  waves  Cplasmons")  to  form  an 
electromagnetic  wave  (“photon")  has  been  studied  theoretically  and 
experimentally  for  a  number  of  years.  (2)  Here  we  briefly  summarize  the 
properties  of  this  three-wave  interaction  as  it  applies  to  the  Two-Stream 
PEM5  concept.  To  recapitulate,  it  is  this  three-wave  interaction  by  which 
some  of  the  steady-state  electron  stream  energy  is  converted  to 
electromagnetic  radiation  at  frequencies  suitable  for  efficient  generation  of 
a  mirror-confined  hot-electron  plasma.  The  DC  power  introduced  through  the 
electron  streams  is  thus  accumulated  and  stored  in  the  anisotropic, 
hot-electon  plasma  in  a  form  from  which  it  can  be  efficiently  transformed 
into  coherent  pulses  of  microwave  radiation. 

As  was  discussed  extensively  by  Harris  (11),  it  is  useful  to  view  the 
three-wave  process  from  a  quantum  mechanical  viewpoint.  Each  wave  can 
then  be  assigned  an  energy,  hoo,  and  momentum,  t>k.  Here  h  is  Planck's 
constant  divided  by  2tt.  Conservation  of  energy  and  momentum  in  the 
three-wave  interaction  thus  requires  that 
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coj  +  c Oy  -  003  And  k  ]  *  k?  -  k3. 

We  have  designated  the  two  electrostatic  waves  with  subscripts  1  and  2, 
and  the  electromagnetic  wave  with  subscript  3.  As  we  have  seen,  the 
fastest  growing  electrostatic  waves  have  frequencies  very  near  the  electron 
plasma  frequency  of  the  background  plasma: 

OOj  =  CO2  cOpe- 

The  wave  vectors  of  the  electrostatic  waves  can  be  nearly  equal  in 
magnitude  but  with  opposite  sign,  so  that  k3  <<  kj  ~  -  k2.  Thus,  the 
electromagnetic  wave  can  be  emitted  with  vanishingly  small  values  of  kM, 
and  will  propagate  radially  outward  away  from  the  spatial  region  occupied  by 
the  streams.  Through  the  three-wave  process  the  electron  streams  serve  as 
a  source  of  bursts  of  microwave  radiation  at  frequencies  near  200,*. 

The  efficiency  with  which  the  DC  electron  stream  energy  is  converted  to 
electromagnetic  power  is  still  somewhat  uncertain.  Several  different 
theoretical  models  have  been  used  to  estimate  this  conversion  efficiency 
(12);  but  the  experimental  measurements  show  substantially  higher 
efficiencies  than  predicted  theoretically.  (2)  The  predicted  parametric 
dependence  of  the  electromagnetic  power  emitted  can  be  seen  from  the 
expression  derived  by  Kaplan  and  Tsytovlch  ( 1 2): 

PM  ~  (2.4TT/\/3)(oOpB2/n0mc3)  J  dk  dJUescop#/ck)2 

OOp^/C 
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Here  is  the  power  density  of  the  electromagnetic  wave,  Ui^  is  the 
saturated  energy  density  of  the  electrostatic  waves  per  unit  interval  of  k, 
and  the  integral  extends  up  to  the  Debye  wave  number,  kD. 

Following  O’Neil  et.  al.  (9)  we  estimate 

~  tWrVn.F3  /Sk. 


where 


5k  ~  k0(nb/f^),/3ir,/2. 

A  rough  estimate  of  the  integral  can  be  obtained  by  assuming  that  has  its 
saturated  value  over  the  interval  (k0  -  $k/2,  k0  +  $k/2)  and  is  zero 
otherwise: 

Pp  ~  0.8tT\/3  (cope2/nemc3)  UJes^Vb/cPtfk. 

The  microwave  power  density  is  thus  predicted  to  vary  roughly  as 
cop*  r\,Eb\/N  (Eb/mc^Kr^/naF  (Vb/cP. 

For  typical  anticipated  system  parameters  the  predicted  power  density  is 
0(1W/cm3]. 
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3.2  EXPERIMENTAL  ARRANGEMENT  AND  OPERATION  OF  THE 

AMPHED  FACILITY 

3.2.1  Description  of  the  AMPHED  facility 

AMPHED  is  a  flexible  CW  magnetic-mirror  facility  comprised  of  six  field 
coils  as  shown  schematically  in  Fig.  7.  The  plasma  confinement  region  is 
enclosed  In  a  cylindrical,  microwave-tight,  aluminum  vacuum  vessel  with 
enough  penetrations  to  accommodate  multiple  microwave  power  sources  and 
plasma  diagnostics. 

The  magnetic  field  configuration  to  be  used  in  the  experiments  proposed 
here  is  produced  by  connecting  the  two  center  coils  in  parallel,  and  then  in 
series  with  the  other  coils  which  are  fed  by  a  200  kVA  DC  power  supply. 
While  the  maximum  overall  mirror  ratio  is  about  2,  the  mirror  ratio  at  the 
cyclotron  resonance  interaction  region  can  be  varied  from  1:1  to  2.1.  The 
location  of  the  ECR  layer  and  the  mirror  ratio  at  this  location  for  a  given 
heating  frequency  can  be  controlled  precisely  by  adjusting  the  D  C.  magnetic 
coil  current.  This  feature  is  particularly  advantageous  for  investigation  and 
optimization  of  electron  cyclotron  heating  mechanisms  which  are  predicted 
to  depend  strongly  on  the  local  magnetic  field  gradient  at  the  location  of 
resonance.  ( 1 4) 

A  plot  of  the  magnetic  field  configuration  is  shown  in  Fig.  8.  (The 
magnetic  field  calculation  has  been  verified  to  an  accuracy  of  0.5%  by 
measuring  the  field  profile  along  the  mirror  axis  using  an  array  of  Hall 
probes.)  There  are  two  noteworthy  features  of  this  magnetic  configuration: 

1.  The  region  of  minimum  field  Is  nearly  uniform  throughout  an 
axial  length  of  80  cm  and  a  radial  extent  of  15  cm,  a  volume  of 
57  liters.  The  average  scale  length  characterizing  the  variation 
of  field  strength  within  this  region  is  greater  than  100  cm. 

2.  Near  the  mirror  throat  the  magnetic  field  strength  increases  to 
twice  its  midplane  value  within  20  cm. 
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The  long  axial  dimension  of  the  AMPHED  facility  also  provides  several 
other  operational  conveniences.  For  example,  the  magnetic  field  near  one  end 
of  the  mirror  configuration  can  be  reduced  locally,  yet  the  field  topology  and 
scale-length  on  the  other  end  can  be  kept  relatively  unaltered.  This 
modification  causes  the  heated  particle  to  gain  additional  parallel  energy 
due  to  the  conservation  of  their  magnetic  moment,  thus  allowing  them  to 
overcome  the  (usually)  positive  ambipolar  potential  well  and  to  be  collected 
by  a  parallel-energy  particle  analyzer,  as  described  in  detail  in  the 
following  section. 

At  present  microwave  power  is  introduced  in  three  different  heating 
configurations,  as  also  shown  in  Fig.  7.  Whistler  waves  are  launched  from  a 
2.45  GHz,  800  watt,  commercial  microwave  oven  coupled  into  the  high  field 
mirror  throat  region  using  a  dielectric-loaded  C-band  waveguide,  terminated 
with  a  tapered  Teflon  slab.  Ordinary  waves,  also  at  2.45  GHz,  are  launched 
from  a  side  port  directly  illuminating  the  flat-field  region  near  the  mirror 
minimum.  The  Teflon  loaded  waveguide  is  oriented  in  such  a  way  that  the 
E-plane  of  the  waveguide  is  parallel  to  the  magnetic  field  lines,  i.e.  E*  II  IT, 
where  T  and!?  are  the  wave  electric  field  and  the  vacuum  magnetic  field, 
respectively.  The  UORH  power,  transmitted  from  a  1 0.1  GHz,  cw  klystron 
amplifier,  is  coupled  to  the  system  through  another  side  port.  This  power  is 
near  the  4th  harmonic  of  the  fundamental  heating  frequency  of  2.45  GHz.  A 
commercial  4-port  hybrid  coupler  is  used  between  the  standard  X-band  wave 
guide  and  the  aluminum  vacuum  vessel  to  improve  the  impedance  matching 
with  the  plasma. 

Plasmas  are  typically  produced  in  argon,  nitrogen  or  hydrogen  gas  at  the 
operating  pressure  of  1CT5  -  lOr3  torr.  Because  significant  plasma  pumping 
generally  occurs  in  this  pressure  range,  the  ionization  gauge  that  monitors 
the  operating  pressure  is  located  in  the  center  region  of  the  vacuum  system 
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to  obtain  more  meaningful  pressure  readings. 

3.2.2  Diagnostics 

The  diagnostics  to  be  used  in  this  investigation  are  indicated  in  Fig.  7. 
The  X-ray  detector  is  a  lead  collimated  commercial  Nal(Tl)  scintillation 
spectrometer  and  Pulse  Height  Analyzer  to  measure  the  hot  electron 
bremsstrahlung  spectrum.  The  detector  is  sensitive  to  X-rays  with  energies 
ranging  from  40  keV  to  the  MeV  range.  The  parallel -energy  particle 
analyzer  is  a  three-grid,  electrostatic,  electron  flux  collector  constructed 
inside  a  1"  diameter  copper  tube  placed  in  the  opposite  mirror  end  from  the 
whistler  wave  launcher.  Because  of  the  small  opening  in  the  separation 
grids,  50  x  50  lines/ inch2  with  80%  transparency,  this  analyzer  is  usually 
only  sensitive  to  the  parallel  energy  of  the  hot-electrons.  Perpendicular 
energy  can  be  measured  by  lowering  the  magnetic  field  strength  near  the 
analyzer  so  that  a  known  fraction  of  the  perpendicular  energy  is  converted 
into  parallel  energy. 

Other  diagnostics,  not  shown  in  Fig.  7,  include  an  axially  moveable  rf 
probe,  a  radially  moveable  Langmuir  probe,  an  azimuthal ly  movable  ionization 
chamber  probe  and  diamagnetic  loops,  etc.  The  rf  probe  signal  is  connected 
to  a  HP8551  spectrum  analyzer  for  detailed  spectral  studies.  The  ionization 
chamber  probe  is  sensitive  to  the  flux  of  electrons  with  energies  greater 
than  1.56  keV.  The  signal  from  this  diagnostic  is  amplified  by  a  sensitive 
electrometer  and  a  synchronous  detector,  and  is  very  useful  in  obtaining  the 
heating  rates  of  the  hot  particles.  Most  of  the  experimental  data  are 
interfaced  to  a  PC/XT  computer  using  board  level  ADC  interface  circuitry. 
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3.2.3  Whistler  wave  injection 

Microwave  power  launched  from  the  high-field  mirror  throat  is  known  to 
excite  whistler  waves  with  the  propagation  vector  F  parallel  to  the 
magnetic  field  lines.  This  was  confirmed  experimentally  by  measuring  the 
spatial  correlation  function,  EAtz]coslzk||(z)],  along  the  field  direction  using 
a  moveable  RF  probe  and  a  microwave  interferometric  technique,  as  shown  in 
Figs.  9a  and  9b.  The  local  wavelength  of  the  observed  wave  decreases  rapidly 
as  it  propagates  toward  the  resonance  layer,  in  qualitative  agreement  with 
the  whistler  wave  dispersion  relation.  (15) 
n2  =  1  -  <*00/(00  -  n  cos  0) 

where  n 2  -  (kc/oo)2  is  the  refractive  index,  a  ~  (cope/oo)2 ,  ft  =  eB/mc,  0  is 
the  angle  between  the  wave  vector  and  the  magnetic  field  lines,  c  is  the 
speed  of  light,  e  is  the  electron  charge,  m  is  the  electron  mass,  oo^  is  the 
electron  plasma  frequency,  and  00  is  the  angular  frequency  of  the  heating 
wave.  We  see  that  as  ft  decreases  and  approaches  00  at  the  ECR  layer,  n 
increases  and  approaches  infinity,  as  observed  in  Fig.  9b.  Notice  that  at  z  ~ 
5  cm,  ft  ~  1.7oo.  Using  the  measured  value  of  a  =  10,  we  have  from  equation 
( 1 ),  n2  =  1  -  I0oo/{oo  -  1.7oo]  ~  15,  or  •  3.8.  This  gives  \  =  K/Z.8  =  3.2 
cm,  in  close  agreement  with  the  measured  wavelength  of  3.3  cm. 
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FIGURE  9a.  Experimental  arrangement  used  for  launching  and  characterizing 
whistler  waves  in  AMPHED. 
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FIGURE  9b:  Experimental  measurement  of  the  whistler-wave  spatial 
waveform,  showing  decreasing  wave  length  as  the  wave 
propagates  toward  the  ECH  layer. 


SECTION  3.3 

Amplification  of  whistler  waves  propagating  through  inhomogeneous, 
anisotropic,  mirror-confined  hot-electron  plasmas 

G.  E.  Guest  and  R.  L.  Miller 
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A  fully  relativistic  local  dispersion  relation  for  whistler  waves  has  been  solved  at  closely  spaced 
points  along  the  magnetic  field  lines  of  a  2:1  magnetic  mirror  in  which  a  highly  anisotropic, 
spatially  inhomogeneous,  hot-electron  plasma  is  confined.  The  limiting  plasma  parameters  for 
convective  (spatial)  growth  have  been  determined  numerically  and  used  to  identify  plasma 
conditions  leading  to  maximum  amplification  of  input  microwave  signals  introduced  in  the 
form  of  whistler  waves.  The  maximum  gain  has  been  evaluated  numerically  for  a  range  of 
values  of  the  hot-electron  plasma  within  which  all  major  stability  criteria  are  satisfied.  Very 


high  gains  ( ~  40  dB)  are  indicated  over  the  entire 

I.  INTRODUCTION 

It  is  well  known  that  whistler  waves  can  be  unstable  in 
mirror-confined  hot-electron  plasmas  because  of  the  pres¬ 
sure  anisotropy  inherent  in  magnetic-mirror  confinement.1 
Indeed,  the  phenomenon  of  unstable  whistler  waves  has  of¬ 
ten  been  invoked  to  interpret  spontaneous  rf  emissions  from 
laboratory2  and  space  plasmas,1  and  has  recently  been  pro¬ 
posed  as  a  way  of  transforming  the  energy  stored  in  mirror- 
confined  hot-electron  plasmas  into  giant  pulses  of  high- 
powered  microwave  radiation.4-'  The  proposed 
transformation  results  from  the  amplification  of  whistler 
waves  propagating  along  magnetic  lines  of  force  and  passing 
through  the  anisotropic,  hot-electron  plasma.  Here  we  in¬ 
vestigate  the  maximum  gain  that  a  whistler  can  undergo  in  a 
spatially  inhomogeneous  plasma,  confined  in  a  magnetic- 
mirror  field  of  finite  spatial  extent. 

Because  the  whistler  instability  mechanism  depends  on 
the  kinetic  properties  of  the  hot-electron  equilibrium,  par¬ 
ticularly  the  distribution  of  electrons  in  velocity  space,  theo¬ 
retical  analyses  of  the  growth  typically  use  idealized,  spatial¬ 
ly  uniform  models  of  the  plasma  in  order  to  treat  the 
velocity-space  issues  realistically.1  Under  these  idealized 
conditions,  growth  of  the  whistler  wave  can  be  either  con¬ 
vective  or  absolute,  depending  on  the  plasma  parameters* 
Absolutely  unstable  or  temporally  growing  waves  are  ex¬ 
pected  to  spread  throughout  the  unstable  plasma  and  grow 
in  amplitude  until  quasilinear  or  nonlinear  processes  balance 
the  destabilizing  conditions.  By  contrast,  the  growth  of  con- 
vectively  unstable  or  spatially  amplifying  waves  may  be  lim¬ 
ited  by  the  finite  transit  time  of  a  wave  packet  across  the 
finite  spatial  extent  of  the  unstable  plasma  as  well  as  by  the 
spatial  nonuniformity  of  the  confining  magnetic  field,  rather 
than  quasilinear  or  nonlinear  mechanisms. 

In  the  present  study  we  use  linearized  theories  to  deter¬ 
mine  the  overall  gain,  G(to),  experienced  by  a  pulse  of  mi¬ 
crowave  power  launched  at  one  end  of  a  mirror-confined, 
hot-electron  plasma  and  propagating  along  the  magnetic 
lines  of  force  in  the  form  of  whistler  waves: 

G(w)  s  |£^,  (ojj  =  +  L)\/\E]n(a>j.  =  -  L) |. 

Here  to  is  the  ( radian )  frequency  of  the  wave  launched  at  the 
axial  position  z  =  -  L  with  rf  electric  field  strength  . 


range  of  beta  investigated 

After  a  single  transit  through  the  plasma,  the  wave  arrives  at 
the  axial  position  z  —  L  with  amplitude  £'ou,  Since  this  pro¬ 
cedure  will  yield  the  desired  gain  only  if  the  unstable  growth 
is  convective,  we  have  carried  out  a  detailed  numerical  eval¬ 
uation  of  the  threshold  conditions  for  absolute  growth  and 
re  .tricted  the  range  of  plasma  parameters  to  lie  entirely  with¬ 
in  the  convectively  unstable  range. 

Apart  from  transients  at  the  beginning  and  end  of  each 
pulse,  the  input  wave  amplitude  is  assumed  to  be  constant  in 
time  throughout  the  pulse;  and  the  spatial  and  temporal  be¬ 
havior  of  the  wave  can  be  determined  using  Fourier  and  La¬ 
place  transformed  equations  in  the  conventional  way  * 
When  inverting  the  transformed  variables  it  is  necessary  to 
displace  the  Laplace  integral  path  to  the  real-ru  axis  while 
deforming  the  Fourier  inversion  path  to  prevent  any  zeros  of 
the  dispersion  relation  from  crossing  the  deformed  path  of 
integration  in  the  complex  k  plane,  where  k  is  the  ( parallel ) 
wavenumber.  As  is  well  known,  this  procedure  is  only  possi¬ 
ble  if  the  whistler  wave  growth  is  convective;  in  the  case  of 
absolute  growth,  two  poles  from  opposite  halves  of  the  com¬ 
plex  k  plane  merge. 

To  ensure  that  we  consider  only  convectively  unstable 
cases,  we  have  evaluated  the  threshold  conditions  for  abso¬ 
lute  (temporal)  growth  for  the  hot-electron  equilibria  under 
investigation.  For  hot-electron  plasmas  created  by  electron 
cyclotron  heating  (ECH),  the  convective/absolute  bound¬ 
ary  can  be  conveniently  represented  in  two  dimensions  by 
fixing  some  of  the  plasma  parameters  in  ways  that  reflect 
fundamental  properties  of  the  ECH  process.  Thus,  for  exam¬ 
ple,  the  total  plasma  density  is  usually  limited  by  microwave 
cutoff*  so  that  (a^/col  <  1.  Here  ru„.  and  <y,,  arc  the  electron 
plasma  and  microwave  frequencies,  respectively.  Moreover, 
the  relative  density  of  hot  electrons  nh/n&b  may  be  limited 
by  the  onset  of  flutelike  instabilities  in  simple  magnetic-mir¬ 
ror  configurations;  although  no  such  limitation  occurs  in 
magnetic  “well’'  configurations  with  favorable  curvature  of 
the  magnetic  lines  of  force.  Even  in  such  a  flute-stable  config¬ 
uration,  however,  the  plasma  pressure  may  be  limited  by 
equilibrium  conditions  such  as  the  so-called  "mirror  mode.” 
and  the  value  of  hot-electron  beta,  P  =  2/z( ,/\/£ may  be 
limited  in  inverse  proportion  to  the  anisotropy.1  Here  B  is 
the  magnetic  intensity  and  //,,  is  the  permeability  of  free 
space. 
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If  no  other  stability  criteria  are  violated,  the  gain  G(  to ) 
can  be  evaluated  up  to  the  threshold  conditions  defined  by 
the  C /A  boundary  in  a  straightforward  way;  this  is  the  main 
objective  of  the  work  described  here.  The  gain  is  then  given 
directly  by 

G(u>)  =  exp^  -  2  J  k,  (coj)dzJ, 

where  k ,  is  the  imaginary  part  of  the  wave  vector  corre¬ 
sponding  to  solutions  of  the  local  dispersion  relation  with 
real  frequency  a>. 

In  order  to  account  for  the  z-dependent  equilibria  con¬ 
fined  in  finite  sized  magnetic-mirror  configurations  we 
adopt  simple  model  equilibria  with  hot-electron  distribution 
functions  mapped  along  magnetic  field  lines  through  adiaba¬ 
tic  invariants.  The  properties  of  these  equilibria  are  summar¬ 
ized  briefly  in  Sec.  II.  A  fully  relativistic  dispersion  relation 
is  then  used  to  evaluate  the  C/A  boundary  for  whistler 
waves.  The  dispersion  relation,  while  fully  relativistic,  can 
give  only  a  local  evaluation  of  the  wave  properties,  since  it  is 
derived  under  the  assumptions  of  spatial  uniformity.  This 
procedure  is  expected  to  yield  a  valid  description  provided 
the  scale  length  for  variation  of  the  plasma  equilibrium  is 
much  greater  than  the  wavelength  of  the  whistler  waves. 
Since  we  are  primarily  interested  in  quasioptical  systems 
with  scale  lengths  much  larger  than  wavelengths,  this  condi¬ 
tion  is  not  particularly  restrictive.  In  effect,  we  are  construct¬ 
ing  a  WKB  solution  from  numerical  solutions  of  the  local 
dispersion  relation  within  plasma  parameter  ranges  for 
which  the  unstable  growth  is  guaranteed  to  be  convective. 
The  numerical  analysis  of  the  dispersion  relation  and  typical 
results  are  recapitulated  in  Sec.  III.  The  gain  is  then  evaluat¬ 
ed  from  real-w  solutions  of  the  dispersion  relation  in  a  way 
that  is  discussed  at  some  length  in  Sec.  IV.  Several  cases  are 
presented  there  in  which  the  temperature  anisotropy  and 
hot-electron  beta  are  varied  so  as  to  move  along  a  particular 
C /A  boundary  (corresponding  to  a  fixed  relative  hot-elec¬ 
tron  concentration,  A  =  30%).  Substantial  values  of  the 
gain  are  indicated  for  hot-electron  plasmas  of  typical  labora¬ 
tory  dimensions. 

The  major  findings  are  summarized  in  Sec.  V,  and  the 
nature  of  these  findings  is  discussed.  In  particular,  the  unre¬ 
solved  issue  of  the  role  of  quasilinear  or  nonlinear  processes 
is  given  some  emphasis  because  of  the  indications  that  the 

_ I 


0,  inside  the  loss  cone. 


If  u,,  is  the  momentum  per  unit  rest  mass  evaluated  at  the 
midplane  z  =  0,  then  in  terms  of  the  constants  of  motion, 

«0||  -  u2  -  =  u\(z)  +  (b-  l)[«?(z)]/£ 

(2) 

and 


maximum  gain  predicted  in  the  linear  theory  \anes  only 
weakly  with  beta  and  remains  large  in  plasmas  with  high 
energy  densities,  where  saturation  levels  of  output  power 
might  be  expected  to  be  very  high. 


II.  THE  MODEL  EQUILIBRIA 

In  this  section  we  summarize  the  way  in  which  the  equi¬ 
librium  properties  of  the  plasma  depend  on  distance  along 
the  magnetic  lines  of  force.  Local  properties  of  electron  cy¬ 
clotron  heated  hot-electron  plasmas  have  been  discussed  ex¬ 
tensively,7  and  a  few  studies  have  attempted  to  model  the 
axial  dependence  in  nonrelativistic  formulations  of  the  hot- 
electron  dynamics.  Here  we  follow  the  conventional  ap¬ 
proach  to  mapping  relativistic-electron  equilibria  along 
magnetic  lines  of  force  assuming  the  invariance  of  the  total 
energy  e,  and  the  magnetic  moment  /j: 

€  =  mc2(y  -  1 )  =  «r  { \1  +  ( trTF  j  -  I  ] 
and 

H  =  (mn])/{2B). 

In  these  expression,  y  is  the  relativistic  Lorentz  factor  and  u 
is  the  momentum  per  unit  rest  mass  with  components  u  and 
t/ii  perpendicular  and  parallel,  respectively,  to  the  magnetic 
field  B.  Here  m  is  the  electron  rest  mass  and  c  is  the  speed  of 
light  in  vacuum 

The  relative  intensity  of  the  steady  magnetic  field  is  as¬ 
sumed  to  vary  along  the  line  of  force  of  interest  as  follows: 

Here  z  =  0  defines  the  so-called  "midplane”  on  which  the 
magnetic  intensity  reaches  its  minimum  value  b  =  1,  while 
z  =  L  defines  the  so-called  "mirror  throat,”  where  B  reaches 
its  maximum  value  b  =  M,  the  mirror  ratio. 

The  cold-electron  component,  which  typically  contains 
the  majority  of  electrons,  will  be  modeled  as  an  isotropic, 
pressureless  fluid  whose  density  is  independent  of  z.  By  con¬ 
trast,  the  hot-electron  component  is  assumed  to  be  aniso¬ 
tropic  and  relativistic  and  described  by  a  distribution  func¬ 
tion  of  the  form’1 


outside  of  the  loss  cone. 


(I) 


r 


=  u]/b. 


(3) 


Here  u0„  and  are  the  components  of  u,,  parallel  and  per¬ 
pendicular,  respectively,  to  the  static  magnetic  field 

A  z-dependent  loss  cone  in  u  space  is  defined  by  setting 
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u\(L)  =  0: 

«S  (*•)  =  =  0 

=  uj  (z)  +  uf  (z)  -  uf  ( z)Bmjtt/B(z ), 

or 


If  we  now  collect  terms  in  u]  (z)  in  the  exponential  factor  we 
have 


where 


a\  (z)  =  (a5„ajjh)/[a4({>  -  1)  +  ajj  ], 
or,  if  tj  s  a /a\, 

a]  (z)  =  T)bal/[r/(b  -  1)  +  1).  (6) 


The  z-dependent  hot-electron  density  is  now  given  by 


n„(z) 


x[*  ^Yej-JLY 
Jiu-VtM-M  \ba]  /  \  a](z)J 


(7) 


and  for  the  simple  bi-Gaussian  distribution,  /  =  0,  we  have 


nh  (z)  =  iryna\  ( z )  or  M  C7  J,u  ’  { 1  +  [ajj/a?(z)] 
X(h/(A/-h)]>',/2 
or 


(8) 


FIG.  I.  Variation  of  relative  hot -electron  density  and  magnetic  intensity,  b, 
with  distance  along  a  magnetic  field  line. 


M*) 

MO) 


oJU)/  (ajj/ain)  (l/(Af  —  D]  V 

a?o  V  1  +  [aj;/a3(z)]  [b/(M~  b)\) 


b  I  (M-b) 

V(b-l)  +  \*\J  (M-\) 


(9) 


An  illustrative  case  is  shown  in  Fig.  1  for  a  mirror  ratio 
M  =2  and  two  different  temperature  anisotropies,  y  =  5 
and  20.  The  z-dependent  relative  magnetic  intensity  b  is 
shown,  together  with  the  z-dependent  relative  hot-electron 
densities  for  the  two  different  temperature  anisotropies.  The 
figure  indicates  the  degree  to  which  higher  temperature  ani¬ 
sotropy  is  associated  with  increased  peaking  of  the  hot-elec¬ 
tron  density  at  the  midplane.  This  property  of  the  equilibria 
will  be  shown  to  have  a  major  impact  on  the  overall  gain  of  a 
whistler  wave  propagating  along  the  magnetic  line  of  force. 
In  quantitative  terms,  the  z  dependence  of  hot-electron  den¬ 
sity  is  given  by 

d  nh(z)  _  db  (  (1-7?)  I  (M-b) 

dz  nh(0)  dz\[\  +  ri(b  -  I)]' yj  (M  -  \) 

— i — ) 

2\I(M  —  1  ){M-b)  1  +  y(b-\) ) 


m  y-y _ '—\ 

dz  V  2(\1  -  1 )  / 


as  b —  1 . 


(10) 


Although  it  is  not  included  in  Fig.  1,  the  temperature  an¬ 
isotropy  also  decreases  with  distance  from  the  midplane  as 

(z)  _  T)b 

crj5  1  +  r j(b  -  1 )  ’ 
so  that 

da]  a\yb'(  1  —  r/) 
dz  [1  +  y(b-  l)]3 

Thus  both  temperature  anisotropy  and  the  relative  hot-elec¬ 
tron  density  decrease  with  distance  away  from  the  midplane, 
•so  that  maximum  amplification  will  occur  at  the  midplane. 

In  contrast  to  this  strong  z  dependence  of  the  hot -elec¬ 
tron  density,  the  cold-electron  density  is  expected  to  remain 
nearly  constant  in  z,  except  for  the  spatial  compression  asso¬ 
ciated  with  the  increasing  magnetic  intensity.  The  real  value 
of  k(w)  is  determined  mainly  by  the  cold-plasma  density 
and  the  magnetic  intensity,  and  is  slowly  varying  along  the 
magnetic  field  lines. 


III.  CONDITIONS  FOR  MAXIMUM  GAIN 

The  maximum  amplification  of  an  externally  launched 
whistler  wave  is  likely  to  be  achieved  in  hot -electron  plasmas 
that  are  just  below  the  threshold  for  absolute  or  temporal 
growth.  This  threshold  is  frequently  referred  to  as  the  C /A 
boundary,  which  is  a  hypersurface  in  a  space  spanned  by  the 
plasma  parameters.  On  one  side  of  this  boundary  only  con¬ 
vective  or  spatial  growth  is  possible,  while  on  the  opposite 
side  of  the  boundary,  absolute  or  temporal  growth  occurs. 
Thus,  in  order  to  explore  the  maximum  gain  that  can  be 
achieved,  it  is  first  necessary  to  define  the  C /A  boundary  in 
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terms  of  equilibrium  parameters  representative  of  physically 
realizable  plasmas.  This  section  summarizes  the  results  of  a 
numerical  determination  of  the  C  /A  boundary  for  unstable 
whistler  waves  in  hot-electron  plasmas  whose  local  proper¬ 
ties  are  modeled  by  the  class  of  infinite,  homogeneous,  rela¬ 
tivistic-electron  distribution  functions  cited  earlier,  namely. 


The  general  methodology  adopted  here  for  determining 
the  C  /A  boundary  has  been  discussed  at  length  by  Briggs'’ 
and  consists  ultimately  in  locating  simultaneous  solutions  of 
the  dispersion  relation  and  its  derivative  with  respect  to  k, 
the  (complex)  wavenumber,  for  real  values  of  the  frequency 
c o .  This  general  methodology  is  implemented  in  our  case 
through  novel  numerical  techniques  appropriate  to  the  class 
of  distribution  functions  used  here.  These  techniques  have 
been  fully  described  by  Guest  and  Miller.9 

It  has  been  well  established  by  many  investigators'  that 
the  growth  rate  of  unstable  whistler  waves  increases  with 
increasing  values  of  the  temperature  anisotropy  and  the  hot- 
electron  beta,  where  beta  is  the  ratio  of  the  kinetic  pressure 
nT  to  the  magnetostatic  pressure  B  72/z„.  For  hot -electron 
plasmas  created  by  electron  cyclotron  heating  (ECH),  sev¬ 
eral  fundamental  aspects  of  the  heating  and  confinement 
processes  limit  the  physically  accessible  values  of  beta  and 
anisotropy  in  ways  that  can  usefully  restrict  the  range  of 
plasma  parameters  to  be  considered  here.7  9  For  example, 
effective  penetration  of  the  ECH  power  typically  limits  the 
maximum  (total)  electron  density  to  values  such  that 
t u <o):M ,  where  ro ,  is  the  angular  frequency  of  the  ECH  pow¬ 
er.  Typically,  the  relative  fraction  of  electrons  heated  to  high 
energies,  A  =  nh/n,  is  much  less  than  unity.  Idealized  theo¬ 
ries  of  flutelike  instabilities  have  been  used  to  estimate  criti¬ 
cal  values  of  A  set  by  thresholds  for  occurrence  of  flutelike  or 
ballooning  modes:  Acnl  —0.3.  Relativistic  effects  appear  to 
increase  this  substantially,  but  we  will  usually  assume 
A<0.5.  In  small  magnetic  mirrors,  nonadiabatic  effects  can 
limit  the  maximum  hot-electron  energy  to  values  such  that 
p<0.03£.  Herep  is  the  (relativistic )  electron  gyroradius  and 
L  is  a  characteristic  scale  length  of  the  magnetic  geometry. 
The  condition  for  equilibrium,  sometimes  referred  to  as  the 
mirror  instability  criterion,  limits  the  product  of  anisotropy 
and  beta,  roughly  as  /?(  rj  —  I)  <  l .  For  the  anisotropic,  hot- 
electron  distribution  functions  used  in  the  present  work,  a 
more  accurate  estimate  of  the  beta  limit  corresponding  to  the 
mirror  criterion  can  be  obtained  from  the  more  general  con¬ 
dition  for  existence  of  equilibria"’: 


,,1  +  i^o 

H0  B  dB 


(12) 


In  the  nonrelativistic  approximation,  p, ,  the  component 
of  plasma  pressure  perpendicular  to  the  static  magnetic  field, 
is  given  by 


pi  =  mn(v\ )  =  mn(v \)  =  mn(v\  /2) 

=  *r I  dut  \  t>,  dViflv^Viimfi. 


(13) 


Hereto,,,  =  v\b  /  (M  —  b)  defines  the  loss  cone  at  the  point  z 
where  the  local  magnetic  intensity  is  B  =  bB„.  Note  that 
y7(n  —0  as  M—  <x .  If  we  use  the  bi-Maxwellian  distribution 
function  corresponding  to  /  =  0K  and  evaluate  the  derivative 
of  p ,  at  the  midplane,  we  obtain  from  the  mirror  criterion 

0  < _ x/A/  —  I  ( 17  - 1  +M  -  \  )'!- _ 

(M-  1 )( l.Sry  1  +M-  1)  —  3r/  78' 

(14) 

where  rjsa\/a\.  We  have  evaluated  the  limiting  value  of 
Pi  ( 77 )  for  3/  as  2  and  displayed  this  boundary  in  Fig.  2.  Note 
that  as  M—  00  this  criterion  approaches 

/?»(*-  1)0. 

In  Fig.  2  we  show  the  value  of  anisotropy  on  the  convec¬ 
tive/absolute  boundary  yc/4  as  a  function  of/7  for  condi¬ 
tions  that  are  appropriate  to  an  ECH-generated  hot-electron 
plasma.  We  have  set  to1  /<o„.  =  1  and  A  =  0.3,  and  varied  y 
to  cover  the  indicated  range  of  values  of  /?  .  The  associated 
values  of  y  are  shown  by  the  curve  labeled  y  —  1 .  Here  y  is 
the  average  of  y  over  the  equilibrium  distribution  function, 
Eq.  ( 5 ),  in  which  only  the  /  =  0  component  is  retained.  The 
approximate  mirror  instability  criteria  /?(?/—  1  )  =  1  and 
t  —  0  are  also  shown  here;  note  that  the  C /A  boundary  lies 
slightly  above  the  equilibrium  limit  for  /?  >  10% -20%.  The 
nonrelativistic  C  /A  boundary  given  by  Iiyoshi  el  at. ' '  is  also 
shown  for  comparison  with  the  fully  relativistic  numerical 
results.  For  a  given  value  of  /?,  the  relativistic  expressions 
indicate  that  slightly  higher  values  of  anisotropy  are  re¬ 
quired  for  absolute  growth. 

The  C/A  boundary  shown  in  Fig.  2  satisfies  the  condi¬ 
tions  for  physical  accessibility  identified  here  for/?  v  10%- 
20%  and  thus  will  be  used  to  estimate  the  maximum  gain 
that  can  be  achieved  at  lower  values  of  beta  The  results  of 


TIG.  2  Variation  of  anisotropy  of  the  C /A  boundary,  y,  ,.  and  average 
hoi -electron  energy.  j\  willi  /?, ,  the  perpendicular  hoi -electron  beta  Here 
idj,  /1  a;',  =  I  and  A  -  0.3.  Also  shown  is  the  mirror  msubiluy  threshold 
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this  evaluation  of  gain  are  described  in  Sec.  IV. 

We  conclude  the  present  section  with  one  further  set  of 
results  on  the  local  wave  properties  of  unstable  whistlers  at 
the  C  /A  boundary  shown  in  Fig.  2.  These  are  displayed  as 
functions  of  / 3 ,  in  Fig.  3,  where  we  see  that  both  the  real  and 
imaginary  components  of  the  propagation  vector,  together 
with  the  wave  frequency,  are  decreasing  functions  of  beta 
This  decrease  reflects  two  underlying  facets  of  the  C/A 
boundary  along  which  the  parameters  are  evaluated,  name¬ 
ly,  that  the  maximum  anisotropy  for  convective  growth  is 
decreasing  with  increasing  values  of  beta,  and  that,  in  the 
present  case,  beta  is  being  increased  by  increasing  the  aver¬ 
age  electron  energy  and  hence  the  relativistic  mass  shift.  The 
relativistic  effects  can  be  removed  by  holding  the  hot-elec¬ 
tron  energy  fixed  and  varying  beta  by  varying  A.  This  proce¬ 
dure  has  been  followed  in  generating  the  results  shown  in 
Figs.  4  and  5. 

Figure  4  again  shows  the  anisotropy  at  the  C  /A  bound¬ 
ary  as  a  function  of  beta,  together  with  the  corresponding 
values  of  A  and  the  mirror  instability  criterion  for  y  =  1.4. 
This  value  of  y  corresponds  to  an  average  hot-electron  ener¬ 
gy  just  above  100  keV  (T  =  111  keV).  Note  that  absolute 
growth  is  only  possible  for  /?i>0.10  corresponding  to 
A  =  0. 1 55,  and  that  the  C  /A  boundary  drops  below  the  mir¬ 
ror  stability  limit  for/?,  >  0. 17  or  A >0.27.  The  wave  param¬ 
eters  along  the  C  /A  boundary  are  shown  in  Fig.  5,  and  here 
reflect  only  the  decreasing  value  of  anisotropy,  as  well  as  the 
increasing  value  of  A.  In  the  range  that  satisfies  the  mirror 
stability  criterion,  the  wave  parameters  decrease  only  slight¬ 
ly,  suggesting  that  relativistic  effects  are  mainly  responsible 
for  the  large  decrease  in  wave  parameters  shown  in  Fig.  3. 

IV.  EVALUATION  OF  MAXIMUM  GAIN 

We  now  estimate  the  maximum  achievable  gain  for  var¬ 
ious  plasma  parameters  by  numerically  evaluating  the  inte- 


*A 


FIG.  3.  Wive  parameters  versus/?,  fort/,  ,^«n  the  C/A  boundary  shown  in 
Fig.  2. 
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FIG.  4.  The  C/A  boundary  for  fixed  y  =  1.4  and  the  relative  hot-electron 
density,  A,  together  with  the  mirror  instability  criterion. 

gral  along  magnetic  field  lines  of  the  imaginary  part  of  the 
propagation  vector: 

In  G(,to)  —  —  2^  k,{io,z)dz.  (15) 

This  procedure  is  valid  provided  the  growth  is  convective. 
We  will  examine  a  sequence  of  cases  in  which  the  plasma 
beta  is  increasing,  and  for  each  value  of  beta  we  will  assume 
the  maximum  anisotropy  for  which  the  whistler  growth  is 
convective.  Since  the  plasma  density  and  anisotropy  are 
greatest  on  the  midplane,  it  is  the  midplane  values  of  plasma 


t± 


Fir.  V  Wave  lur  ndns  along  llu-  C/A  boundary  for  the  cases  shown  in 

Fig.  4. 
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and  wave  parameters  on  the  C  /A  boundary  that  determine 
the  maximum  anisotropy  and  the  corresponding  wave  fre¬ 
quency  for  each  value  of  beta.  The  anisotropy  and  density 
are  readily  determined  at  each  point  along  the  magnetic  field 
line;  and  the  dispersion  relation  is  then  solved  locally  to  ob¬ 
tain  the  complex  propagation  vector  at  that  point  for  the 
frequency  specified  at  the  midplane.  The  gain  integral  is  then 
evaluated  numerically. 

We  have  carried  out  this  evaluation  of  the  gain  for 
=5%,  10%,  15%,  20%,  and  25%  in  the  case  shown  in  Fig. 
2  =  1,  A  =  0.30,  and/ =  0.  The  magnetic  geometry 

assumed  was  a  fixed  M  —  2  magnetic  mirror;  and  thus  no 
attempt  was  made  to  relate  the  equilibrium  value  of  the  an¬ 
isotropy  to  the  shape  of  the  magnetic  field.  We  return  to  this 
point  in  Sec.  V.  The  midplane  parameters  for  the  cases  exam¬ 
ined  are  tabulated  in  Table  I.  The  series  was  terminated  at  P, 
=  25%  since  the  mirror  instability  criterion  was  not  satis¬ 
fied  at  the  C  / A  boundary  for  greater  values  of  /?, . 

The  last  column  of  Table  I  gives  the  value  of 


The  gain  is  then  given  by 

G(<u)  =  exp[4 (Lu>ctU/C)I  ] . 

For  typical  conditions,  the  applied  microwave  power  reso¬ 
nates  with  the  second  harmonic  of  the  local  electron  gyro  fre¬ 
quency  at  the  midplane.  Thus 

<w«o  = 

and  we  can  estimate  the  scale  factor  as  follows: 

ALiO'M  4irLf„  _  4 irL 

c  c  A„ 

where  A,,  is  the  free-space  wavelength  of  the  microwave 
heating  power.  In  order  of  magnitude,  10)  for 

typical  existing  experimental  devices.  The  indicated  values 
of  gain  are  then,  again,  in  order  of  magnitude, 

G(«)~0(  10*)~40dB. 

In  addition  to  the  peak  growth  rates  shown  in  Table  I  we 
can  readily  estimate  the  effective  bandwidth  by  calculating 
the  total  gain  for  whistler  waves  of  varying  frequencies  prop¬ 
agating  through  a  particular  anisotropic,  mirror-confined, 
hot-electron  equilibrium.  For  definitiveness,  we  again 
choose  the  plasma  parameters  to  be  near  the  convective/ 


FIG.  6.  The  frequency-dependent  gam  factor,  •  j ( A.  t Ulz/l..  for 
whistler  waves  of  varying  frequency  in  hoi -electron  equilibria  mill  tem¬ 
perature  anisotropies  near  the  conveciise/absolute  boundary  and/?  V", . 
KWi .  and  20%.  I’lasma  and  wave  parameters  for  these  three  cases  arc  listed 

in  Table  I. 


absolute  boundary  for  the  center  frequency  of  the  unstable 
frequency  band.  Illustrative  cases  of  /?,  =  5%,  10%,  and 
20%  are  shown  in  Fig.  6,  where  the  dimensionless  gain  fac¬ 
tor. 


is  plotted  against  the  frequency  normalized  to  the  center  fre¬ 
quency.  Note  that  the  bandwidth  increases  as/?  increases, 
i.e.,  as  the  average  hot-electron  energy  increases.  The  full- 
width  at  half-maximum  varies  from  roughly  0.15  at  /? 
=  5%  to  more  than  0.25  at  /?,  =  20%. 

Thus  the  linear  theory  suggests  that  very  high  gain  is 
possible.  It  should  be  emphasized  that  the  linear  theory  can- 


TABLE  I.  Summary  of  the  wave  and  plasma  parameters  at  the  C/A  boundary  for  varying  values  of/?,  For  all  these  cases,  -  I  and  i-0  JO 


0,1%) 

V<  /A 

*i 

-  k,c/b>.,„ 

y 

r  k  c  dz 

i  /. 

5 

15.26 

0.6998 

I  623 

0.7267 

1.093 

0.076 

10 

11.83 

0  5769 

1.316 

0.5964 

1.197 

0073 

15 

963 

04840 

1  093 

0.5150 

1  3125 

0.070 

20 

7.972 

04095 

09357 

0.4487 

1.4375 

0067 

25 

6.698 

0  3484 

0.810*1 

0.3942 

l  572 

0  064 

30 

5.724 

0.2977 

0.7097 

0.3485 

1.717 

35 

4.988 

02562 

0.6280 

0  3106 

1.870 

40 

4,426 

0.2220 

0  5617 

0.2787 

2.031 
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not  provide  any  indication  of  the  maximum  output  power  or, 
equivalently,  the  saturation  value  of  the  input  power.  These 
must  be  obtained  from  quasilinear  or  nonlinear  models  of  the 
energy  transformation  process. 

It  is,  in  fact,  most  striking  that  the  gain  decreases  only 
very  slowly  with  increasing  beta;  one  would  expect  that  the 
maximum  output  power  would  be  proportional  to  beta.  As 
we  remarked  earlier,  the  fact  that  the  spatial  growth  rate 
decreases  with  beta  is  in  part  attributable  to  the  somewhat 
arbitrary  assumption  of  fixed  A.  Beta  is  then  proportional  to 
the  average  energy  of  the  hot  electrons,  and  relativistic  ef¬ 
fects  reduce  the  growth  rates.  Perhaps  more  important  to  the 
weak  beta  dependence  of  the  gain  is  the  fact  that  at  low  beta 
the  C/A  boundary  permits  very  high  anisotropy.  High  an¬ 
isotropy,  in  turn,  leads  to  very  large  spatial  growth  rates  near 
the  midplane.  Although  lower  anisotropy  does  lead  to  less 
peaking  of  the  density  near  the  midplane,  the  growth  is 
greatly  reduced  away  from  the  midplane  by  the  increased 
local  magnetic  field  strength:  o>/u>cc  is  too  small  to  permit 
large  local  growth  rates. 

V.  CONCLUSIONS 

In  Sec.  IV  we  demonstrated  that  the  amplification  of  an 
externally  launched  whistler  wave  can  exceed  40  dB  over  a 
relatively  wide  range  of  values  of  beta,  5%  </?,  <25%,  for 
typical  ECH  plasma  parameters:  w^/col,  =  1,  A  =  nh/n 
=  0.3,  and  1 . 1  <  y<  1 .6.  The  corresponding  limiting  values  of 
temperature  anisotropy  were  15>?;r/<4  >6.  These  values  of 
anisotropy  are  much  larger  than  expected  in  a  typical  ECH 
plasma  confined  in  a  2-.1  magnetic  mirror.  Dandl4  has  pro¬ 
posed  the  use  of  adiabatic  compression  to  increase  rj  to  the 
C/A  boundary  value.  We  have  not  attempted  to  include 
adiabatic  compression  in  the  present  studies.  The  gain  calcu¬ 
lation  included  the  strong  spatial  variation  in  hot-electron 
density  and  anisotropy,  together  with  the  spatial  variation  in 
the  underlying  2: 1  magnetic-mirror  field.  We  conclude  from 
these  findings  that  strong  amplification  is  possible  for  repre¬ 
sentative  hot -electron  equilibria  in  the  linear  regime;  i.e.,  for 
input  signals  below  some  power  level  at  which  quasilinear  or 
nonlinear  effects  will  saturate  the  growth  of  the  whistler 
wave.  Clearly  the  linear  theory  cannot  provide  an  estimate  of 
this  critical  parameter;  some  type  of  quasilinear  calculation 


of  the  saturation  of  whistler  growth  is  required  to  address 
this  issue. 

Existing  quasilinear  studies  of  unstable  whistler  waves 
in  infinite  homogeneous  plasmas  have  indicated  that,  within 
that  highly  idealized  model,  roughly  10%  of  the  hot-electron 
energy  can  be  transformed  into  wave  energy  before  the 
growth  rate  is  reduced  to  zero.  If  these  estimates  provide  a 
valid  indication  for  finite,  inhomogeneous,  mirror-confined 
plasmas,  it  would  appear  entirely  feasible,  by  operating  at 
higher  values  of  beta  in  high  field  strength  magnetic  fields  to 
achieve  very  high  output  power  levels.  In  fact,  it  may  well  be 
that  the  infinite  homogeneous  quasilinear  model  provides  an 
underestimate  of  the  saturation  power,  since  resonant  elec¬ 
trons  can  escape  a  real  magnetic  mirror  after  giving  a  sub¬ 
stantial  fraction  of  their  perpendicular  kinetic  energy  to  the 
wave.  In  effect,  a  nearly  constant  level  of  anisotropy  is  auto¬ 
matically  maintained  in  an  actual  magnetic  mirror,  whereas 
the  infinite,  homogeneous  quasilinear  model  allows  the  ani¬ 
sotropy  to  decrease  as  a  result  of  the  wave  grow  th. 
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ABSTRACT.  A  one-dimensional,  steady-state,  relativistic  Fokker-Planck  model  of  electron  cyclotron  heating 
(ECH)  is  used  to  analyse  the  heating  kinetics  underlying  the  formation  of  the  two-component  hot-electron  plasmas 
characteristic  of  ECH  in  magnetic  mirror  configurations.  The  mode)  is  first  applied  to  the  well  diagnosed  plasmas 
obtained  in  SM-1  and  is  then  used  to  simulate  the  effective  generation  of  relativistic  electrons  by  upper  off-resonant 
heating  (UORH),  as  demonstrated  empirically  in  ELMO.  The  characteristics  of  unstable  whistler  modes  and  cyclotron 
maser  modes  are  then  determined  for  two-component  hot -electron  plasmas  sustained  by  UORH.  Cyclotron  maser 
modes  are  shown  to  be  strongly  suppressed  by  the  colder  background  electron  species,  while  the  growth  rates  of 
whistler  modes  are  reduced  by  relativistic  effects  to  levels  that  may  render  them  unobservable,  provided  the 
hot-electron  pressure  anisotropy  is  below  an  energy  dependent  threshold. 


1.  INTRODUCTION 

Early  experiments  on  electron  cyclotron  heating 
(ECH)  of  plasmas  confined  in  magnetic  mirrors  1 1  ] 
using  CW  microwave  power  at  a  single  frequency 
created  an  unusual  plasma  containing  two  distinct 
populations  of  electrons.  The  majority  of  the  electrons 
remained  at  a  temperature  that  was  typically  below 
100  eV,  while  a  small  fraction  was  heated  to  average 
energies  above  100  keV.  Rather  broad  regimes  of 
stable  operation  were  found  to  be  limited  by  the  onset 
of  several  dominant  modes  of  instability,  such  as  the 
mirror  mode,  the  flute  mode  and  unstable  whistler 
waves.  Later  experiments  in  the  ELMO  facility  [2] 
showed  that  the  energy  stored  in  the  hot-electron 
component  could  be  greatly  increased  by  using  upper 
off-resonant  heating  (UORH),  in  which  microwave 
power  at  a  frequency  above  the  cold-electron  gyro- 
frequency  was  added  to  the  fundamental  ECH  power. 
In  addition  to  increasing  the  relative  density  of  the  hot 
electrons,  UORH  also  enhanced  the  stability  of  the 
plasma,  even  for  values  of  beta  in  excess  of  50%  (beta 
is  the  conventional  ratio  of  plasma  kinetic  pressure  to 
magnetostatic  pressure). 

The  stable,  high-beta,  hot-electron  annulus  revealed 
by  the  ELMO  experiments  led  to  the  successful  use  of 
these  plasmas  ('ELMO  rings’)  to  stabilize  a  scalar 
pressure  core  plasma  confined  in  the  ELMO  Bumpy 


Torus  device  (see,  for  example.  Ref.  [3]):  this  stimu¬ 
lated  widespread  interest  in  the  fundamental  properties 
of  these  remarkable  plasmas  (see,  for  example. 

Ref.  (4)).  Recently  proposed  applications  for  these  hot- 
electron  plasmas  include  stabilization  of  localized 
modes  in  heliacs  [5]  and  tokamaks  [6],  as  well  as 
energy  storage  for  subsequent  conversion  to  high- 
power  microwave  pulses. 

It  has  required  the  efforts  of  many  workers  to 
evolve  reasonably  satisfactory  theoretical  models  of  the 
steady  state  equilibrium  produced  in  ECH  experiments. 
The  basic  characteristics  of  the  interaction  of  individual 
electrons  with  microwave  electric  fields  were  derived 
very  early  (for  a  discussion  of  early  theoretical  work, 
see  Ref.  [7])  and  useful  formulas  were  obtained  for  the 
resulting  heating  rate  under  a  wide  range  of  conditions 
(for  a  description  of  later  work,  see  Ref.  [81).  Sub¬ 
sequent  development  of  geometrical  optics  codes  (see, 
for  example.  Ref.  [9J)  made  it  possible  to  describe  the 
propagation  and  absorption  of  microwaves  in  specified 
plasma  media.  Progress  was  accelerated  by  the  formula¬ 
tion  of  a  relativistic  Fokker-Planck  theory  of  ECH  by 
Bernstein  and  Baxter  [10].  Hamasaki  and  co-workers 
(Il|  took  moments  of  the  Fokker-Planck  equation, 
using  trial  distribution  functions,  in  an  attempt  to 
predict  the  power  required  for  formation  of  the 
ELMO  rings.  Other  workers  [12]  developed  numerical 
techniques  for  a  direct  solution  of  the  Fokker-Planck 
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equation;  these  are  now  implemented  in  large  computer 
codes  [13].  Ultimately,  these  codes  are  expected  to 
provide  a  detailed,  self-consistent  picture  of  high-beta 
hot-electron  plasmas  formed  by  ECH. 

Nonetheless,  there  is  an  immediate  need  for  a 
reduced  model  that  can  be  used  to  search  over  the 
wide  range  of  parameters  characterizing  an  ECH 
experiment.  It  is  not  practical  to  carry  out  suci 
search  with  two-dimensional  Fokker-Planck  cooes; 
this  would  require  broad,  time  consuming  parameter 
searches  to  overcome  the  general  lack  of  experimental 
measurements  of  such  essential  parameters  as  the 
strength  of  each  RF  electric  field  component  present  in 
the  plasma.  We  have  addressed  this  need  by  evolving  a 
Fokker-Planck  model  for  the  steady  state  distribution 
of  electrons  in  energy;  this  model  includes  in  a  heuristic 
way  the  phenomena  that  depend  essentially  on  the 
electron  orientation  in  a  gyrotropic,  two-dimensional 
velocity  space.  This  model  has  been  described  in 
considerable  detail  in  Ref.  [14], 

In  order  to  model  a  complex  physical  situation  in  a 
relatively  simple  way,  the  power  balance  issues  affecting 
the  formation  of  a  relativistic  electron  tail  have  been 
considered  primarily  and  not  the  particle  balance 
issues,  except  in  so  far  as  a  significant  associated 
energy  loss  occurs,  as  in  the  case  of  non-adiabatic 
scattering  of  high  energy  electrons  into  the  loss  cone. 
For  a  more  detailed  justification  of  this  approach,  see 
the  Appendix.  Thus,  we  have  chosen  to  conserve 
particles  while  accounting  for  the  energy  loss  by 
converting  the  non-adiabatic  loss  term  into  a  dynamical 
friction  term,  as  described  in  Appendix  C  of  Ref.  [14], 
This  choice  has  the  desired  physical  effect  of  sharply 
reducing  the  hot-electron  distribution  function  for 
energies  above  the  threshold  for  violation  of  adiabatic 
invariance,  as  is  clearly  evident  in  Figs  la  and  lb.  The 
hot-electron  distribution  in  energy  is  shown  in  Fig.  la 
for  plasmas  typical  of  those  reported  in  Ref.  [15]. 
Figure  lb  shows  the  four  rates  accounted  for  in  the 
model;  the  quasi-linear  heating  rate,  the  Coulomb 
thermalization  rate,  the  synchrotron  radiation  cooling 
rate  and  the  energy  loss  rate  due  to  non-adiabatic 
scattering. 

Here,  we  first  apply  this  reduced  model  to  the  SM-1 
experiments  [15]  in  order  to  demonstrate  the  success  of 
this  model  in  describing  the  hot-electron  plasmas 
studied  in  SM- 1 .  The  results  of  this  demonstration  are 
presented  and  discussed  in  Section  2.  They  clarify  a 
number  of  features  of  ECH  experiments  that  had  not 
been  fully  understood  previously.  In  particular,  the 
role  of  microwave  electric  fields  parallel  to  the  static 
magnetic  field  (ordinary  mode  of  propagation)  is 


FIG  la.  Electron  distribution  in  energy  for  a  case  typical  of  the 
single-frequency  heating  experiments  reported  in  Ref  {151. 


FIG.  lb.  Characteristic  rates  entering  into  the  heating  kinetics 
(from  top):  electron  heating  rale  due  to  ECH  power;  electron 
thermalization  by  Coulomb  collisions  with  colder  electrons: 
synchrotron  radiation;  and  energy  loss  by  non-adiabatic  scattering 
of  hot  electrons  into  the  loss  cone. 


shown  to  be  critical  to  an  effective  formation  of  the 
hot-electron  component.  Furthermore,  by  exploiting 
the  sensitivity  of  the  hot-electron  heating  rate  to  the 
parallel  microwave  electric  fields,  we  are  able  to  use 
the  reduced  Fokkcr- Planck  ECH  model  to  demonstrate 
the  impact  of  non-adiabatic  scattering  [16],  which  has 
been  conjectured  for  some  time  to  limit  the  maximum 
temperature  of  the  hot  electrons  in  most  ECH  experi¬ 
ments  to  date. 
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We  then  apply  the  present  model  to  the  UORH 
experiments  in  ELMO  [2],  where  the  model  displays 
the  same  striking  increases  in  stored  energy  with 
increasing  UORH  power  that  were  observed  in  ELMO. 
Finally,  the  model  is  used  to  obtain  high-beta  hot- 
electron  equilibria  suitable  for  subsequent  stability 
studies.  Although  the  reduced  model  cannot  provide 
information  on  pressure  anisotropies  or  other  two- 
dimensional  properties  of  the  equilibria,  it  does  relate 
the  relative  temperatures  and  densities  of  the  two 
electron  species  and  provides  a  much  larger  degree  of 
internal  consistency  than  was  heretofore  available. 

In  Section  3,  we  analyse  the  stability  of  the  hot- 
electron  plasmas  to  transverse  waves  propagating  along 
the  static  magnetic  field,  namely  whistler  waves  [17] 
and  the  cyclotron  maser  modes  (see,  for  example, 

Ref.  [18a]  and  the  excellent  set  of  references  cited 
therein;  for  more  recent  results,  see  Ref.  [18b]).  It  has 
been  known  for  many  years  [19]  that  relativistic  effects 
generally  reduce  the  growth  rates  of  unstable  whistler 
waves,  but  specific  results  depend  on  the  assumed 
equilibrium.  Here  we  use  the  Fokker-Planck  model  to 
obtain  equilibria  that  are  consistent  with  ECH  kinetics 
and  demonstrate  the  relativistic  suppression  of  whistler 
growth  rates.  The  cyclotron  maser  modes  are  found  to 
be  completely  stabilized  by  the  cold-electron 
component. 

Finally,  we  note  in  passing  that  our  understanding 
of  flute  modes  in  high-beta  hot-electron  plasmas  has 
not  been  considered  here  and  remains  incomplete, 
although  fully  electromagnetic  kinetic  models  of  line 
tying  [20]  have  helped  clarify  finite-length  effects. 
There  is  an  active  effort  to  predict  the  stability  of  these 
high-beta  plasmas  to  flute-like  modes  [21]  which  we 
hope  will  be  aided  by  more  detailed  descriptions  of 
realistic  equilibria.  Discussion  and  conclusions  are 
found  in  Section  4. 


2.  FORMATION  OF 

HOT-ELECTRON  PLASMAS  USING  ECH 

In  this  section  we  describe  the  results  of  applying 
the  reduced  Fokker-Planck  ECH  model  to  ECH 
experiments  in  the  mirror  devices  SM-1  [15]  and 
ELMO  [2]  before  using  the  model  to  obtain  high-beta 
equilibria  for  stability  analysis.  The  Fokker-Planck 
model  is  described  in  great  detail  in  Ref.  [14]  and 
references  cited  there. 

Briefly,  we  use  a  one-dimensional,  steady-state, 
bounce-averaged  Fokker-Planck  model  of  the  ECH 
process.  Since  our  primary  concern  is  with  the  thermali- 


TABLE  1.  SM-1  REFERENCE  CASE  (Ref  [  15]) 


Stored  energy 
Radius  of  annulus 
Radial  width 
Axial  length 
Annulus  volume 
Mirror  ratio  on  axis 
Coil  separation 
Mirror  ratio  at  r0 
Mirror  ratio  at  resonance 
ECH  power 
Microwave  frequency 

Background  plasma 
density 
temperature 

Hot-electron  component 
density 
temperature 


W  =  20  J 
r,(  =  10  -i  cm 
2.3  =  5.8  cm 
L'2  =  12  anu  18  cm 
Vol  =  12  x  I0’  cm’ 
M„  =  2.2 
L  =  7 1  cm 
M  =  3 

m„,  =  : 

r>  «=  1  k\V 
f  «  9.5  GH/ 

n  =  ( I  -  5 1  x  10 1  ‘  cm 
T  »=  10-50  eV 

n„  =  3  x  10 !l  cm 
T„  »  350-400  ke\ 


zation  of  a  relatively  small  concentration  of  energetic 
electrons,  we  adopt  a  linearized,  test-particle  frame¬ 
work.  The  relativistic  Choi')  electrons  arc  then 
regarded  as  a  distinct  species  whose  spatial  density  n„ 
is  much  smaller  than  the  densities  of  the  other  plasma 
species.  All  of  these  ‘field’  species  are  assumed  to  be 
isotropic,  thermal,  and  at  a  much  lower  temperature 
than  the  hot  electrons.  This  model  determines  the 
response  of  a  small  relative-density,  high-energy 
electron  tail  to  the  incident  microwave  power  b> 
including  the  dominant  dynamical  processes:  quasi- 
linear  heating,  thermalization  by  Coulomb  scattering  on 
cooler  electrons,  synchrotron  radiation,  and  direct  loss 
by  non-adiabatic  scattering.  Specifically,  the  model 
predicts  the  electron  distribution  in  energy  resulting 
from  RF  electric  fields  with  specified  values  of 
frequency,  amplitude  and  polarization,  incident  upon  a 
background  plasma  of  specified  density  and  tempera¬ 
ture.  confined  in  a  prescribed  mirror-like  magnetic 
field. 

We  take  as  our  SM-1  reference  case  the  20  J 
equilibrium  discussed  at  length  in  Ref.  [15]  and 
summarized  in  Table  I.  As  mentioned  earlier,  there  are 
no  direct  measurements  of  the  strength  of  the  micro- 
wave  electric  fields,  and  we  must  scan  over  all 
plausible  values  of  both  components.  E,  and  E,  We  do 
this  by  specifying  the  ratio  E|/E„  of  the  O-mode  and 
X-mode  field  strengths  and  searching  out  the  values  of 
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FIG.  2a.  Perpendicular  component  of  the  microwave  electric  field, 
,  average  energy  of  the  hot  electrons,  Th,  and  total  RF  power 
density,  Pu.  versus  background  plasma  temperature,  T,,  for  an 
assumed  ratio  of  O-mode  to  X-mode  field  strengths,  Et/E  L  =  2, 
in  the  SM-1  configuration  summarized  in  Table  I.  The  relative 
hot-electron  density  is  njn  =  0  1. 


FIG.  2b.  Variation  of  E  L,  Th  and  RF  power  density  Pu  with 
E,/Em  for  a  background  plasma  temperature  Tr  =  30  eV  and 
nk/n  m  01. 


Ex  required  to  sustain  the  observed  relative  hot- 
electron  density,  nh/n  =  0.1;  this  yields  the  observed 
hot-electron  average  energy  Th  and  the  RF  power 
density  required  to  sustain  the  equilibrium.  Since 
the  microwave  power  is  launched  primarily  in  ordinary 
modes  and  is  multiply  reflected  from  the  cavity  walls, 
we  assume  an  essentially  isotropic  distribution  of  wave 


vectors  k  for  the  ECH  waves.  An  illustrative  case  is 
shown  in  Fig.  2a,  which  displays  EA,  Th  and  PRF 
versus  the  background  electron  temperature  Tc  for 
E^Ej.  =  2.  The  characteristic  scale  length  of  the  static 
magnetic  field  has  been  set  at  L|  =  30  cm  for  this 
case,  a  choice  that  we  will  discuss  at  length  in  the 
following.  Note  that  the  RF  power  density  required  to 
sustain  this  equilibrium  is  around  20  mW-cm"3, 
accounting  for  roughly  250  W  of  RF  power.  This  does 
not  include  the  fundamental  ECH  power  necessary  to 
sustain  the  background  plasma.  If  the  energy  expended 
in  creating  each  ion-electron  pair  is  roughly  50  eV,  we 
estimate  this  fundamental  ECH  power  requirement  at 
25  mW'Cm'5  for  an  ambient  hydrogen  gas  pressure  of 
5  x  10  "6  torr.  Since  the  cold-plasma  volume  is 
estimated  at  25  L  in  Ref.  [15],  as  much  as  625  W  of 
fundamental  ECH  power  may  be  required  simply  to 
sustain  the  background  plasma  and  only  250  W  may  be 
required  to  sustain  the  hot-electron  component.  The 
variation  of  Ex,  Th  and  PRF  with  the  ratio  E,/Ex  is 
shown  in  Fig.  2b  for  a  background  electron  tempera¬ 
ture  Tc  =  30  eV.  Note  the  rapid  rise  in  RF  power 
density  for  E|/EA  >  10,  in  which  PRF  increases  five 
times  while  the  average  energy  of  the  hot  electrons 
remains  almost  constant.  The  saturation  of  T„  despite  a 
substantial  increase  in  RF  power  indicates  the  rapid 
increase  in  the  non-adiabatic  scattering  rate.  If  these 
calculations  are  repeated  for  different  values  of  the 
magnetic  scale  length  L|,  we  obtain  the  striking  results 
shown  in  Fig.  3a.  It  is  especially  instructive  to  display 
the  L|  dependence  of  the  saturated  value  of  the  hot- 
electron  average  energy  and  the  corresponding  values 
of  the  electron  gyroradius: 


=  1.7  X  10’  Gem  V«(2  +  e)/B 

where  e  =  Thill/moC:  and  B  is  the  magnetic  intensity. 
Since  we  require  p  <  kL|  for  adiabatic  confinement  of 
electrons,  we  can  evaluate  k  by  plotting  pM1  versus  L(. 
The  value  of  k  obtained  from  Fig.  3b  is  0.06.  which 
agrees  reasonably  well  with  the  conventional  result  cited 
by  Uckan  [16],  From  the  values  obtained  for  Th  u,  it 
appears  that  L(  =  24  cm  is  a  reasonable  value  of  the 
magnetic  scale  length  in  the  hot-electron  region  of 
SM-1.  From  a  parabolic  fit  to  the  vacuum  magnetic 
intensity  we  would  estimate  L*  =21  cm,  in  reasonable 
agreement.  Note  that  this  is  much  less  than  Rc.  the 
radius  of  curvature  of  the  magnetic  field  lines  in  this 
region,  which  is  estimated  to  be  around  60  cm. 

For  large  values  of  L»,  particularly  at  higher 
magnetic  fields,  synchrotron  radiation  becomes  a 
dominant  loss  mechanism  in  the  model.  This  loss  is 
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FIG.  3a.  Th  versus  for  different  values  of  the  magnetic  scale 
length  L\. 


FIG.  3b.  Saturated  value  of  Th  m  and  corresponding  value  of  the 
average  hot-electron  gyroradius  pm  versus  L,.  The  solid  line  is  for 
pm  m  0.06  L\. 


circumvented  in  the  case  discussed  here  by  specifying 
values  of  L|  similar  to  the  ELMO  (high  field)  experi¬ 
ments:  L|  =  10  cm  (ELMO). 

Another  important  feature  of  the  dependence  of  T„ 
on  Ei/E^  is  the  rapid  initial  increse  in  Th,  with  PRF 
constant  or  even  decreasing.  This  behaviour  makes 
clear  the  substantial  benefit  of  increased  microwave 
integrity  of  the  vacuum  chamber.  If  E|/EA  &  10  can 
be  sustained  in  a  suitably  designed  cavity  where  wall 
reflections  do  not  rapidly  convert  ordinary  modes 
to  extraordinary  modes,  Th  can  be  increased  by 
roughly  100  keV  or  more.  Since  the  value  of  E| 
required  to  achieve  this  level  of  heating  is  approxi¬ 
mately  10J  V-cm'1,  while  the  power  density  dissi¬ 
pated  in  the  plasma  is  around  1  Wem~\  the  effective 
quality  factor  must  exceed  300.  The  ELMO  device  [2] 
utilized  a  very  high-Q  cavity  together  with  UORH 
power  to  obtain  exceptionally  efficient  heating  of 
relativistic  electrons.  The  frequency  of  the  UORH 
power  is  chosen  to  be  above  the  cold  electron  gyro- 
frequency  throughout  the  entire  cavity.  This  power 
is  therefore  absorbed  primarily  by  the  relativistic 
electrons  and  only  negligibly  by  the  background 
electrons.  Even  for  the  hot  electrons,  the  absorption 
rate  is  small:  thus,  if  the  cavity  has  a  large  quality 
factor,  the  electric  field  strength  can  increase  to  large 
enough  values  to  render  the  heating  rate  (and  hence  the 
absorption  rate)  appreciable.  Since  heating  rates  are 
quadratic  in  electric  field  strength,  they  can  be  large 
enough  to  sustain  high  energy  electron  populations  at 
levels  below  100  V-cm'1. 

The  Fokker-Planck  simulation  of  early  ELMO 
experiments  on  UORH  [2]  are  she  a  n  in  Fig.  4.  Here 
the  fundamental  ECH  is  at  a  frequency  of  10.6  GHz 
while  the  UORH  source  has  a  frequency  of  35.7  GHz. 
These  results  compare  favourably  with  the  empirical 
results  displayed  in  Fig.  4  of  Ref.  [2].  The  background 
plasma  density  and  temperature  are  2  x  10"  cm'5 
and  50  eV,  respectively,  and  we  have  arbitrarily  set 
E|/Ex  =  10  for  both  sources.  In  the  absence  of 
UORH  power,  the  relative  density  of  hot  electrons  is 
roughly  2.5%,  sustained  by  an  RF  power  density  of 
less  than  6  mW-cm'5.  This  relative  density  increases 
to  64%  and  the  hot-electron  beta  approaches  50%  as 
the  UORH  power  is  increased. 

Having  calibrated  our  model  on  SM-1  and  ELMO, 
it  is  of  interest  to  see  how  much  hot-e1 .  ctron  energy 
can  be  stored  in  comparably  sized  devices  by  optimizing 
the  parameters.  By  using  UORH  at  higher  frequencies 
(60/90  GHz)  in  a  device  that  is  similar  in  physical 
scale  to  ELMO,  we  can  obtain  Fokker-Planck 
equilibria  with  average  hot-electron  energies  above 
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FIG.  4.  Average  hol-electron  energy,  Th,  relative  hot-electron  den¬ 
sity,  njn,  and  plasma  beta,  0,  versus  power  density  from  the 
UORH  source,  Poorh-  The  parameters  describing  the  background 
plasma  and  the  heating  configuration  are  given  in  the  text. 


moCJ.  The  relativistic  effects  on  plasma  stability  can  be 
clearly  examined  in  this  hot-electron  temperature 
range.  In  these  studies,  we  have  again  specified  the 
background  plasma  density  at  n  =  1  x  10 13  cm'3  and 
varied  the  background  electron  temperature  Tc  from 
250  to  3000  eV.  We  assume  that  optimal  background 
electron  temperatures  can  be  maintained  by  a  suitable 
level  of  fundamental  ECH  power,  depending  mainly  on 
the  confinement  properties  of  the  configuration.  For 
the  Fokker-Planck  model  of  the  hot-electron  population 
we  reduce  the  fundamental  ECH  electric  field  strength 
to  a  negligible  level:  10  ~3  times  the  UORH  field 
strength.  We  then  vary  Tc  over  its  range  to  determine 
the  values  of  UORH  power  density  PRf  required  to 
sustain  specified  relative  densities  of  hot  electrons.  The 
Fokker-Planck  model  evaluates  Th  for  specified  values 
of  the  various  system  parameters,  such  as  the  magnetic 
field  scale  length  L|,  here  set  at  10  cm,  as  in  ELMO. 
These  results  are  displayed  in  Fig.  5.  Note  that,  since 
these  curves  are  deduced  for  constant  hot-electron 
densities  Ah  =  nh/n,  the  resulting  value  of  beta  is 
simply  porportional  to  Th.  The  RF  power  density 
exhibits  a  broad  minimum  that  moves  from  500  eV  to 
800  eV  as  the  hot-electron  density  ratio  is  increased 
from  10%  to  40%.  The  associated  beta  values 
approach  50%  and  the  hot-electron  average  energy  is 
around  550  keV.  These  are  the  equilibria  whose 
stability  will  be  studied  in  the  following  sections. 


3.  STABILITY  ANALYSIS 
FOR  TRANSVERSE  WAVES 
PROPAGATING  PARALLEL  TO 
THE  STATIC  MAGNETIC  FIELD 

As  discussed  in  Section  2,  the  first  task  addressed  in 
this  work  was  to  obtain  high-beta  equilibria  from  the 
reduced  one-dimensional  Fokker-Planck  model  of  the 
ECH  process.  Guided  by  the  empirical  results  of 
Dandl  et  al.  [2]  on  UORH  in  the  ELMO  experiments, 
we  found  a  corresponding  regime  in  the  Fokker-Planck 
analysis  that  yields  high-beta,  relativistic-electron 
equilibria  ideally  suited  to  a  stability  analysis,  namely 
beta  values  around  50%  and  relativistic-electron 
temperatures  above  500  keV.  Our  one-dimensional 
Fokker-Planck  model  cannot  provide  any  information 
on  the  temperature  anisotropy  that  is  characteristic  of 
hot-electron  plasmas  formed  by  ECH  and  which  plays 
a  dominant  role  in  the  stability  of  these  plasmas  |22] 
Fully  two-dimensional  Fokker-Planck  ECH  codes 
are  now  available  and  will  eventually  provide  two- 
dimensional  equilibria.  These  numerical  equilibria  can 
be  expressed  analytically  as  a  superposition  of  the 
equilibrium  distribution  functions  described  by  Guest 
and  Dory  [23]: 
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FIG  5.  RF  power  density  required  and  beta  values  achieved  m 
maintaining  specified  values  of  relative  hot-electron  density  versus 
background  electron  temperature  Tr.  The  background  plasma 
density  is  I0"  cm ' 1  and  the  ECH  frequencies  are  60  and  90  GHz. 
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fn(u:,  ui) 


n 

T3/2aia|f 


f,(r,  u,  t)  = 


dt' 


x  exp 


x 


kui  \  d(<j_  +  _kui  dfp 
yw  /  <3uf  yw  3ui 


where  f  =*  0,  1,  2,  ....  and  u  =  p/m  is  the  momen¬ 
tum  per  unit  rest  mass.  Here,  we  examine  our  one¬ 
dimensional  equilibria  for  stability  against  transverse 
waves  propagating  parallel  to  the  static  magnetic  field, 
with  the  temperature  anisotropy  ij  *  Tx/T,  and  the 
loss-cone  index  l  set  at  arbitrary  values.  We  single  out 
the  whistler  mode  [17,  19,  24]  and  the  cyclotron  maser 
mode  [18,  24]  because  of  the  relative  stringency  of  the 
stability  criteria  for  electrostatic  and  whistler  modes,  as 
discussed  by  Guest  and  Sigmar  [22],  and  because  the 
cyclotron  maser  mode  was  omitted  from  consideration 
in  Ref.  [22]. 

In  what  follows,  we  describe  an  efficient  and 
very  rapid  numerical  technique  for  solving  the  fully 
relativistic  dispersion  relation  for  right-hand,  circularly 
polarized,  transverse  waves  propagating  parallel  to  the 
equilibrium  magnetic  field.  The  derivation  of  the 
dispersion  relation  is  sketched  here  for  completeness ; 
more  detailed  discussions  are  readily  available  [25]. 

The  relativistic  Vlasov  equation  for  electrons  is 

(df/dl)  +  (u/y)(df/<3r) 

-  (e/m)[E  +  (u  x  B)/y](3f/du)  =  0 


x  u^fE^e'10  +  E.e'°) 

where  <t>  is  the  polar  angle.  In  the  uniform  magnetic 
field  assumed  here,  the  unperturbed  electron  orbits  are 
given  by  z  =  Zo  +  U|t/y  and  <t>  =  +  eBt/7m,  with 

U|  and  uA  constant.  Choosing  Zo  =  0  at  t  =  0.  we 
obtain 

ft(z  =  0,  u,  t  =  0)  =  -iV2ux(e/m) 


/  ku,  N 

\  dfo 

ku, 

dfo  ' 

'  3u; 

4-  - 

a: 

dui 

/  E,e  10  E-e'° 

x(  +  - 

\kU|  -  7u>  -  eB/m  kui  -  yw  +  eB/m 

From  Faraday’s  and  Ampere's  laws. 

V  X  (V  X  E)  =  -Hofdld t)[j  +  e0OE/dt)] 
or 

k2E  -  (k-E)k  —  (wJ/c2)E  =  iw/xj 


where  y2  =  1  +  u2/c2  is  the  Lorentz  factor  and  the 
velocity  is  v  =  u/y.  The  perturbed  distribution 
function,  associated  with  a  wave  whose  electric  and 
magnetic  fields  are  E  and  B  =  k  x  E/oi,  is  given  by 


For  purely  transverse  waves  propagating  parallel 
to  the  equilibrium  magnetic  field,  with  components 
Et  *  (E,  ±  iE,)/V2,  this  expression  becomes 


Since  k-E  =  0  for  transverse  waves,  E  and  j  must 
satisfy 

l(k2c2/w2)  -  1]E  =  (i/oKo)]* 

For  the  right-hand,  circularly  polarized  component,  the 
resulting  dispersion  relation  is 


X  - 

kui  -  yw  +  eB/m 


NUCLEAR  FUSION.  Vot  21,  No  )  (I9U) 


50 


425 


GUEST  and  MILLER 


For  the  class  of  equilibrium  distribution  functions 
discussed  by  Guest  and  Dory  [23],  the  numerator  of 
the  integrand  can  be  rewritten  conveniently  as 


+  fot  dfo 
70)  3uj 


Since  the  Lorentz  factor  7  appeals  in  the  resonant 
denominator  and  since  we  are  concerned  with  strongly 
relativistic  plasmas  with  7  2  2,  it  is  necessary  to 
evaluate  the  integral  numerically.  We  now  consider 
integrals  of  the  general  form 


I(«.k)  =  f"duie-u*,/ai  f”  du,  e-u'/“' 

Jo  J  -• 

x  F(U|,uf,M,k) 
ku,  —  70)  4-  eB/m 


For  complex  values  of  w  and  k,  there  may  be  zero, 
one  or  two  solutions  to  the  relativistic  resonance 
condition  ku(  -  70)  +  eB/m,  which  we  denote 
by  ui,  (i  =  1,  2).  We  denote  each  solution  for 
U|(u,k,uJ)  by  w  and  the  corresponding  values  of  7 
by  7?  —  v  1  +  w2/c2  +  u2/c2.  Now  consider  the 
following  integral: 


r. 


du  F(w„u2  ,o),k) 

*  kU|  -  7,u  +  eB/m 


a,F(w 
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dy 
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where  Z  is  the  conventional  plasma  dispersion  function 
[26] .  From  this,  we  obtain: 


I(u,  k)  = 


du-  e'u-/a- 


V7T 
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7,u*  -  eB'm 
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J  0 


F(U|,  ui,  a),  k) 
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F(w„  u~ ,  w,  k)  \ 
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With  the  singularities  removed  in  this  manner,  the  u, 
integral  can  be  evaluated  using  the  Gauss-IIermite 
quadrature  [27]: 


dui 


e'u'/o'  3(ui)  =  w«...  2(ui.,) 


Similarly,  the  ux  integral  can  be  evaluated  using  the 
Gauss-Lag uerre  quadrature  [26] 

[”du2  e'u'/Q‘>(ui)  =  £  WL  l*(u;.,) 

J  i>  1 

To  evaluate  the  dispersion  relation  for  given  values  of 
u)  and  k,  we  then  choose  a  value  of  s  (s  =  12  typically) 
and  for  each  of  the  s  prescribed  values  of  u;, 

(i  =  1 . s)  we  evaluate  the  corresponding  root  or 

roots  of  the  relativistic  resonance  condition,  u,  ,  =  w„ 
such  that 


krw,  -  7,o),  +  eB/m  =  0 


The  Gauss-Hermite  formula  is  then  used  to  calculate 
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FIG.  6.  Density  dependence  of  cyclotron  maser  mode  growth 
rate  for  fixed  values  of  beta.  The  plasma  parameters  are  specified 
in  the  text. 


A  standard  routine  is  used  to  evaluate  the  plasma 
dispersion  function,  and  the  Gauss-Laguerre  quadrature 
is  used  to  evaluate  the  u?  integral: 


like/loss-cone  aspect  of  the  mirror-confined  relativistic- 
electron  equilibria  [28].  Of  particular  interest  for  the 
present  studies  is  the  stabilizing  effect  of  a  cold- 
electron  component  on  the  cyclotron  maser  mode. 

Since  our  Fokker-Planck  ECH  equilibria  consist  of  a 
minority  population  of  relativistic  electrons  in  a 
majority  population  of  much  colder  background 
electrons,  cyclotron  maser  modes  are  anticipated  to  be 
strongly  suppressed.  This  is  illustrated  in  the  results 
shown  in  Fig.  6,  where  we  fix  the  beta  value  at  50% 
and  100%,  and  the  relative  density  of  hot  electrons 
at  30%.  The  hot-electron  distribution  has  a  loss 
cone  index  (  =  2,  and  a  temperature  anisotropy 
r)  *  Tx/T|  =  3.  The  peak  temporal  growth  rate  w, 
normalized  by  the  (non-relativistic)  electron  gyro- 
frequency  o)«  is  plotted  against  the  total  electron 
density,  specified  through  the  ratio  of  electron  plasma 
and  gyrofrequencies,  (u)pr/uct)2.  For  the  parameters  of 
interest,  the  peak  growth  rates  are  less  than  10" 3  times 
the  gyrofrequency.  Figures  7a  and  7b  demonstrate  the 
fundamental  dependence  of  the  cyclotron  maser  modes 
on  the  loss-cone  index  f,  rather  than  the  temperature 
anisotropy.  For  these  cases,  the  hot-electron  beta  is 
fixed  at  30%,  the  total  density  at  =  w*c/2.  and  the 
relative  density  of  hot  electrons  at  50%. 
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FIG.  7a.  Demonstration  of  the  independence  of  the  cy  clotron 
maser  mode  growth  rate  with  respect  to  temperature  anisotropy. 


With  this  algorithm  for  evaluating  the  velocity-space 
integrals,  we  can  use  standard  root-finding  routines  to 
solve  the  dispersion  relation  for  fully  relativistic  hot- 
electron  plasmas.  In  this  way  we  can  evaluate  the 
threshold  conditions  and  growth  rates  for  unstable 
whistler  and  cyclotron  maser  modes  in  the  high-beta 
hot-electron  plasmas  obtained  from  our  Fokker-Planck 
model.  As  has  been  pointed  out  by  Chu  and  Hirshfield 
[18a],  whistler  waves  with  phase  velocities  below  the 
speed  of  light  are  destabilized  by  the  temperature 
anisotropy,  while  the  cyclotron  maser  mode, 
corresponding  to  waves  with  phase  velocities  greater 
than  the  speed  of  light,  are  destabilized  by  the  beam- 


FIG.  7b.  Dependence  of  the  cyclotron  maser  mode  growth  rate  on 
the  toss-cone  index  t. 
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FIG  8.  Stabilization  of  the  cyclotron  maser  mode  by  background 
cold  electrons.  The  relative  density  of  hot  electrons  is  A,  =  njn. 


Finally,  Fig.  8  shows  how  the  growth  is  suppressed 
by  decreasing  the  relative  density  of  hot  electrons,  Ah. 
In  this  particular  case,  stabilization  is  complete  for 
Ah  S  40%;  here,  0  =  30%  and  =  0.5.  Tsang 

[24]  considered  this  effect,  but  only  for  large  values  of 
Ah  where  he  did  not  see  complete  stabilization. 

By  comparison,  the  whistler  waves  can  have  growth 
rates  that  are  larger  by  two  orders  of  magnitude  if  the 
anisotropy  exceeds  a  critical  value,  which  is  about 
three  for  parameters  characteristic  of  hot-electron 
plasmas  formed  by  ECH.  Indeed,  whistler  growth  rates 
can  be  reduced  to  low  values  mainly  by  increasing  the 
average  energy  of  the  hot  electrons  and  decreasing 
the  temperature  anisotropy.  These  two  effects  are 
illustrated  in  Figs  9a  and  9b,  in  a  single-species  hot- 
electron  plasma  with  (  -  0,  jj  =  3  and  0  =  30%. 

The  average  perpendicular  energy  is  expressed  as 
=  Tfc/mc2;  in  Fig.  9b,  =  1.  Finally,  we  evaluate 

the  whistler  growth  rates  in  a  two-component 
equilibrium  obtained  from  the  Fokker-Planck  ECH 
model:  T„  =  527  keV,  Ah  =  30%,  Tc  =  800  eV, 

=  0-5  and  0  =  31%.  The  temperature 
anisotropy  has  been  varied  (i?  =  3,  6,  10,  20)  to  show 
the  rapid  decrease  in  growth  rate  for  decreasing 
anisotropy,  as  displayed  in  Fig.  10a.  For  the  weakly 
unstable  whistler  waves  the  growth  is  convective.  It  is 
thus  of  interest  to  evaluate  the  spatial  amplification 
rates  k„  as  illustrated  in  Fig.  10b. 


4.  DISCUSSION  AND  CONCLUSIONS 

We  have  demonstrated  several  points  in  this  paper. 
First  of  all,  using  the  classical  Fokker-Planck  model  of 
ECH,  it  is  possible  to  reproduce  all  of  the  features  of 
the  SM-1  experimental  results  that  can  be  compared 
with  the  reduced  one-dimensional  model.  In  particular, 
there  is  no  need  to  invoke  collective  processes,  such  as 


FIG.  9a  Whistler  nunle  growth  rates  for  fixed  antsotrop\ .  rj  - 
and  different  values  if  the  average  hoi -electron  energy  7 he  plasma 
parameters  are  specified  in  the  text. 


FIG.  9h  Whistler  nunle  growth  rates  for  fixed  average  hot- 
electron  energy,  t  i  -  I,  and  different  wines  of  the  temperature 
anisotropy 
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unstable  whistler  waves,  to  provide  a  self-consistent 
picture  of  the  hot-electron  distribution  in  energy  and 
the  RF  power  density  required  to  sustain  it.  This 
property  of  the  present  Fokker-Planck  model  is  in 
contrast  to  the  moment  approach  of  Hamasaki  et  al. 
[11],  as  one  might  expect  from  the  non-Maxwellian 
character  of  the  hot-electron  equilibrium  distribution. 
An  interesting  possibility  is  suggested  by  these  results, 
namely  that  a  different  class  of  test  distribution  func¬ 
tions,  similar  to  those  obtained  here,  could  be  used  in 
the  moment  approach  to  obtain  useful  analytic  descrip¬ 
tions  of  ECH  equilibria. 

The  numerical  description  obtained  in  the  present 
work  also  contributes  to  clarifying  the  role  of  non- 
adiabatic  scattering  in  limiting  the  hot-electron  energy 
and  it  reveals  in  quantitative  terms  the  importance  of 
the  microwave  integrity  of  the  cavity  in  which  the 
ECH  is  carried  out.  For  single-frequency  heating  it  is 
apparently  necessary  to  have  a  high  enough  cavity 
quality  factor  to  permit  the  ordinary  mode  field 
strength  to  exceed  by  a  large  factor  the  extraordinary 
mode  electric  field  strength.  Similarly,  UORH  can 
only  be  expected  to  yield  the  high  energy  densities 
observed  in  ELMO  if  the  UORH  field  strengths  can 
approach  values  of  the  order  of  102  V-cm'1,  while  the 
corresponding  absorbed  power  densities  remain  around 
1  W -cm'5  or  less. 

These  results  are  of  practical  importance  in  designing 
experimental  ECH  approaches  to  achieving  relativistic- 
electron  plasmas;  but  they  also  suggest  that  the 
reduced  Fokker-Planck  model  can  provide  a  useful 
simplified  description  of  other  plasmas  in  which  ECH 
may  occur,  such  as  magnetospheric  plasmas. 

For  our  immediate  purposes,  the  reduced  Fokker- 
Planck  model  provides  a  detailed,  self-consistent 
description  of  hot-electron  plasmas  that  can  be  analysed 
for  stability.  In  keeping  with  a  large  body  of  empirical 
observations,  we  find  that  whistler  waves  are  likely  to 
be  unstable  unless  the  temperature  anisotropy  is 
reduced  and  the  hot-electron  energy  is  increased. 
Cyclotron  maser  modes  are  shown  to  be  suppressed  by 
the  cold  background  electrons  that  normally  make  up 
the  majority  species  in  an  ECH  plasma. 

In  a  recent  PhD  thesis,  Gamer  [29]  has  studied  the 
characteristics  of  microinstabilities  in  an  ECH  plasma 
produced  in  a  magnetic  well  configuration  by  single- 
frequency  heating.  Gamer  finds  that  the  hot-electron 
component,  whose  temperature  is  400  keV,  is  remark¬ 
ably  stable;  the  observed  unstable  whistler  waves  are 
apparently  driven  by  a  distinct  group  of  warm  electrons 
at  a  temperature  around  2  keV.  The  third  electron 
component  is  a  cold  population  at  a  temperature 


FIG  10a.  Whistler  mode  temporal  growth  rates  for  a  Fokker- 
Planck  ECH  equilibrium  with  ad  hoc  values  of  the  temperature 
anisotropy. 


FIG.  10b.  Whistler  mode  spatial  growth  rates  and  propagation 
vector  for  the  Fokker-Planck  ECH  equilibrium  with  if  =  J. 

around  100  eV.  The  cold-,  warm-  and  hot-electron 
densities  are  roughly  equal  and  are  in  the  range  of 
(1-2)  x  10"  cm-3.  Garner’s  thesis  contains  an 
interesting  speculation  regarding  the  observed  stability 
of  the  hot-electron  group  that  is  not  predicted  by  his 
analyses,  namely  that,  if  stochastic  heating  can  occur 
at  energies  above  the  predicted  threshold  for  super- 
adiabatic  electron  behaviour,  the  resulting  distribution 
of  hot  electrons  might  indeed  be  stable  at  high  temper¬ 
atures.  Our  kinetic  model  of  the  relativistic  heating 
process  assumes  that  stochastic  heating  is  maintained 
by  the  overlap  of  relativistic,  Doppler-shifted 
resonances,  and  our  analysis  of  the  resulting  equilibria 
demonstrates  the  relativistic  stabilization  of  whistler 
waves  as  well  as  the  cold-electron  stabilization  of  the 
fast-wave  branch  of  the  electromagnetic  mode.  The 
distinct  warm-electron  component  can  be  minimized 
or  even  eliminated  by  the  use  of  O-mode  heating  or 
UORH,  as  shown  by  the  present  Fokker-Planck 
results. 
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APPENDIX 

In  brief,  our  picture  of  particle  balance  in  mirror- 
confined  ECH  plasmas  is  as  follows: 


A-l.  Low-energy  processes 

At  energies  e  5  100  eV,  particle  balance  is  main¬ 
tained  by  three  dominant  processes,  namely  (i)  ioniza¬ 
tion  of  ambient  neutral  gas,  (ii)  free  flow  of  cold  ions 
along  the  magnetic  lines  of  force  and  out  the  ends  of 
the  magnetic  mirror,  and  (iii)  electrostatic  confinement 
of  electrons  in  an  ambipolar  potential  well  whose 
depth  $  is  governed  by  the  electron  heating  rate  i  and 
the  cold-ion  escape  rate. 

Equating  the  ion  and  electron  confinement  times 
yields 


For  typical  values  of  the  ECH  electric  field, 

Erf  «  10-20  Vcm"\  the  central  magnetic  field, 

B0  =  3  kG,  the  electron  density,  n*  *  5  x  10"  cm'3, 
the  overall  mirror  ratio,  M  =  2,  and  the  mirror  ratio 
at  the  test  electron  turning  point,  M,  =  1.7,  the  criti¬ 
cal  electron  energy  for  runaway,  e,,  can  be  just  above 
100  eV.  (Here,  K  is  the  elliptic  integral  associated 
with  the  electron  bounce  motion.) 

Simultaneously,  with  ECH,  small-angle  Coulomb 
collisions  scatter  heated  electrons  into  the  loss  cone  at 
a  rate  given  by 

_6  n*  In  A  eV3'2  cm3 


The  corresponding  rate  of  pitch-angle  diffusion 
associated  with  ECH  is  given  by 


(ir)  * 15  * 10 


e$  =  Le/Vj 


where  v,  is  the  cold-ion  flow  velocity  and  L  is  the 
length  of  the  plasma  column.  Furthermore,  the  escape 
of  these  electrons  and  cold  ions  is  balanced  by  ioniza¬ 
tion,  so  that 

Vj/L  —  nQ(<rKmiaiK)n^e) 


where  no  is  the  density  of  neutral  gas  in  the  plasma. 
Under  typical  conditions,  the  characteristic  rates  for 
these  three  processes  are  around  5  x  104  s'1. 

A-2.  Intermediate-energy  processess 

In  the  energy  range  100  eV  £  «  £  10  keV,  the 
dominant  dynamical  processes  are  Coulomb  scattering 
and  ECH.  The  Coulomb  scattering  causes  both 
slowing-down  and  pitch-angle  scattering  of  heated 
electrons.  For  energy  balance,  the  slowing-down  must 
be  offset  by  ECH  for  energies  above  a  relatively  low 
critical  energy  t  X  »  3e4>: 

I<Rf(*|)I  ~  l<Coul(tl)l 

Using  rudimentary  formulas  for  each  rate,  we  obtain 


where  Mrci  is  the  mirror  ratio  at  the  resonant  surface. 
Since  0RF  *  « -l,  while  0Coo,  =  e~3  2,  RF-driven  pitch- 
angle  diffusion  may  dominate  at  high  energies,  t  > 

*  10  keV,  typically.  In  the  intermediate  energy  range, 

« ,  £  t  <  tt.  Coulomb  scattering  is  typically  the 
dominant  loss  process.  The  resulting  pitch-angle 
displacement  of  an  electron  with  initial  pitch  angle  0, 
can  be  estimated  from 

A0  =  f'2  ^2  df 

J,t  «RF 

Since  eRF  is  almost  independent  of  energy  in  this 
range,  we  have 

A0  _  1.5  x  10~6  n.  In  A 

8,  8}  <RF 


«Ek  (  M-l  V”  ,  (  M,-l  \ 
2BoM,  {  M,-l  /  \  M-l  ) 


n.  In  A 


tWin  cm3 


1.3  x  105  t!a 


eV3/2  cm3 


f<2 

dt 

»  M 


-3/2 


3  x  10  6  nc  In  A 

8*  f*RF 


(ff,:  ~  ff':)  - 


eV>:  cm’ 
s 
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Typically,  (A0/0,)  <  1 ,  provided  the  heating  rate  is 
sufficiently  large,  for  example  eRF  <  107  eV-s'1. 

Some  of  the  runaway  electrons  will  indeed  be  lost, 
especially  for  £  2  6,,  but  the  fractional  hot-electron 
density  obtained  from  energy  balance  alone  provides  a 
useful  estimate  of  an  upper  bound  on  this  ratio.  Note 
also  that  in  the  one-dimensional  Fokker-Planck  model 
we  define  the  hot-electron  density  as  consisting  only  of 
electrons  with  energies  above  a  specified  value,  which 
is  usually  chosen  as  20  keV.  For  «  >  20  keV,  the  RF- 
driven  pitch-angle  diffusion  typically  exceeds  pitch- 
angle  scattering  due  to  Coulomb  collisions;  the  loss 
rate  of  heated  electrons  in  this  energy  range  is  esti¬ 
mated  to  be  very  small.  Since  we  specify  the  total 
electron  density  in  the  one-dimensional  model  and 
restrict  the  hot-electron  component  to  those  electrons 
for  which  Coulomb  scattering  is  negligible,  we  expect 
the  relative  hot-electron  density  obtained  from  energy 
balance  considerations  to  be  a  valid  upper  bound  and  a 
close  estimate. 

A-3.  High-energy  processes 

In  the  energy  range  t  >  10  keV,  Coulomb  scattering 
pitch-angle  diffusion  into  the  loss  cone  can  typically  be 
offset  by  RF-driven  pitch-angle  diffusion  in  which 
heated  electrons  diffuse  towards  the  resonant  surface 
(where  uRF  *  NeB/ym,  N  =  1,2,  3,...).  Heating  at 
the  second  and  higher  harmonics  of  the  relativistic- 
electron  gyrofrequency  plays  a  major  role  in  the  main¬ 
tenance  of  efficient  stochastic  heating  and  causes 
further  pitch-angle  diffusion  towards  the  midplane  of 
the  magnetic  mirror,  0  =  *72.  Particle  losses  in  this 
energy  region  are  very  small,  and  energy  balance  is 
dominated  by  relativistic  ECH  and  synchrotron  radia¬ 
tion  up  to  the  energy  at  which  electrons  cease  to  be 
adiabatic.  The  escape  of  a  very  small  number  of  very 
high  energy  electrons  has  little  effect  on  particle 
balance,  but  can  completely  offset  continued  heating  to 
higher  energies. 

In  summary,  we  suggest  that  the  justification  of 
separating  energy  balance  and  panicle  balance  condi¬ 
tions  rests  on  the  smallness  of  (A0/0),  as  estimated 
above,  as  well  as  on  the  unusual  nature  of  the  RF- 
driven  pitch-angle  diffusion  that  tends  to  drive  heated 
electrons  back  towards  the  resonant  surface.  We  have 
attempted  to  minimize  the  impact  of  panicle  balance 
processes  on  our  one-dimensional  Fokker-Planck 
model  by  regarding  only  electrons  with  energies  above 
20  keV  as  ‘hot  electrons’  and  specifying  the  total 
density  ab  initio. 
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4.  INTERACTIONS 

Progress  In  achieving  the  Phase  1  objectives  was  enhanced  through 
discussions  with  Dr.  Joseph  Santoru,  of  the  Hughes  Research  Laboratories,  and 
Dr.  Marlene  Rosenberg,  of  Jaycor.  Dr.  Santoru  is  active  in  an  experimental 
program  to  characterize  and  exploit  three-wave  interactions  in  a  novel 
microwave  source  concept;  Dr.  Rosenberg  has  provided  a  major  portion  of  the 
theoretical  support  for  the  Hughes  project.  Both  scientists  provided 
pre-publication  reports  describing  detailed  experimental  and  theoretical 
results  regarding  the  characteristics  of  the  three-wave  interaction  and 
various  parasitic  phenomena.  Their  work,  supported  by  AFOSR,  has  been 
valuable  in  demonstrating,  at  least  in  unmagnetized  plasmas,  the 
effectiveness  of  the  three-wave  interaction  in  generating  electromagnetic 
radiation  that  is  ideally  suited  for  creating  the  magnetic-mirror  confined 
hot-electron  plasmas  needed  for  amplification  of  whistlers. 

5.  PATENT  DISCLOSURE 

A  patent  application  relating  to  this  concept,  "Method  and  Apparatus  for 
Producing  Intense  Microwave  Pulses"  Serial  No.  07/293,307,  was  filed  on 
January  4,  1 989. 
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